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ABSTRACT 

A  case  of  explosive  extratropical  cyclogenesis  in  the 
North  Pacific  is  analyzed  employing  guasi-Lagrangian  diag- 

j 

nostic  techniques  id  isobaric  coordinates.  The  First  GARP 
Global  Experiment  (FGGE)  Level  Ill-b  data  set  from  the  Euro¬ 
pean  Canter  For  •  Medium-ra  nge  Weather  Forecasting  (ECMWF)  is 
used  in  a  synoptic  investigation  and  mass  and  vorticity 
budget  diagnostic  evaluation  of  storm  development. 

Bid-tropospheric  positive  vorticity  advection  (PVA)  acts 
in  concert  with  low-level  thermal  advection  to  enhance 
surface  layer  organization  during  the  initial  periods. 
Explosive  development  occurs  under  upper-level  zonal  flow 
when  the  surface  system,  developing  along  the  low-level 
baroclinic  zone,  moves  under  the  front  left  quadrant  of  a 
strong  jet  streak.  Stability  in  the  lower  troposphere 
decreases  steadily  as  an  intense  two-layer  mass  circulation 
is  established.  The  divergence  term  provides  the  largest 
contribution  to  the  vorticity  tendency  during  explosive 
development  due  to  substantial  surface  layer  convergence  and 
upper-level  divergence. 
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I .  INTROD  OCTIOH 


Extratropical  maritime  cylones  restrict  commercial, 
industrial  and  military  use  of  the  sea.  The  more  intense 
storms  which  form  over  very  small  time  periods  are  espe¬ 
cially  dangerous  to  navigation  due  to  the  effects  of  the 
associated  wind,  waves  and  visibility  restrictions.  Accu¬ 
rate  and  reliable  prediction  of  these  atmospheric  processes 
has  keen  a  goal  for  many  decades.  During  the  past  few  dec¬ 
ades,  significant  advances  in  forecasting  these  synoptic- 
scale  atmospheric  circulation  features  have  been  realized 
through  the  development  of  numerical  weather  prediction 
(NWP)  models. 

While  current  numerical  models  provide  satisfactory 
accuracy  in  forecasting  most  synoptic-scale  processes,  they 
perform  poorly  in  predicting  the  predominantly  cold-season 
explosive  genesis  of  maritime  cyclones.  Sanders  and  Gyakum 
(1980)  define  explosive  cyclogenesis  as  an  extratropical 
cyclone  which  undergoes  central  pressure  falls  of  at  least  1 
mb/h  for  24  h.  Accurate  forecasting  of  these  meteorological 
events  is  limited  by  the  scarcity  of  observations  over 
oceanic  areas,  which  inhibits  the  model  initialization  pro¬ 
cess,  and  by  the  smaller  spatial  and  temporal  scales  of  the 
explosively  developing  cyclone. 
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Accurate  forecasting  of  these  maxitine  events  is  of 
great  significance  to  commercial,  industrial  and  military 
organizations  which  routinely  use  the  oceans  in  the  trans¬ 
action  of  business  and  in  the  fulfillment  of  commitments. 
Several  recently  publicized  examples,  specifically,  the  less 
of  the  Russian  trawler  Metros  toy  (Mariners  weather  Leg, 

1979)  ,  and  the  sinking  of  the  oil  rig  Ocean  Ranger  (Lemoyne, 
1982),  illustrate  only  a  minute  fraction  of  the  total  loss 
of  life  and  operational  hazards  associated  with  these 
intense  metacrolcgica 1  events. 

Explosive  deepening  is  a  characteristic  of  the  vast 
majority  of  the  most  intense  cyclones  (Sanders  and  Gyakum, 

1980) .  In  connection  with  associated  severe  weather  and  the 
failure  of  numerical  forecast  models  to  accurately  predict 
these  events,  it  is  of  paramount  importance  that  the 
processes  which  contribute  to  the  formation  of  these  vigor¬ 
ous  storms  be  examined  and  fully  understood.  Several  recent 
studies  have  focussed  attention  on  detailed  analyses  of 
these  processes.  Roman  (1981)  completed  a  mass  budget  study 
and  conant  (1982)  completed  a  circulation  and  angular  momen¬ 
tum  budget  study  of  a  case  of  east-coast  cyclogenesis. 
Intercompariscns  cf  these  results  with  conclusions  derived 
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in  this  thesis  provide  valuable  insight  into  the  formative 
mechanisms  of  maritime  explosive  cyclogenesis  and  lead  tc  a 
general  quantitative  understanding  of  the  forecast  problem 
associated  with  these  intense  systems. 

This  thesis  is  part  cf  a  broadly-based  research  effort 
attempting  tc  understand  better  the  formative  processes 
inherent  in  the  genesis  of  maritime  cyclones.  Specifically, 
this  study  concentrates  on  a  clearly  defined  structured 
approach  to  the  analysis  of  a  case  of  explosive  cyclogenesis 
over  the  Ncrth  Pacific  Ocean  during  the  period  12-17  January 
1979.  The  mass  and  circulation  (vorticity)  budgets  are 
examined  using  quasi-Lagrangian  diagnostics  ( QLD }  as  devel¬ 
oped  by  Jchnscn  and  Downey  (1975  a  and  b)  in  order  to  deter¬ 
mine  the  role  cf  dynamic  forcing  in  the  budget  volume.  The 
local  contxibuticn  of  vorticity  sources  and  sinks  and  the 
vertical  and  horizontal  mass  and  vorticity  transports  are 
specifically  explored. 

Idantif icaticn  of  these  processes  is  uniquely  interest¬ 
ing  in  that  the  cyclone  under  study  consists  of  the  classic 
frontal  wave  development  from  strong  low  level  baroclinic- 
ity.  The  distinguishing  feature  is  the  absence  of  a  strong 
migratory  short  wave  aloft,  which  categorizes  this  storm  as 
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a  Petterssen  type-A  svstem,  as  opposed  to  the  type-B  fea¬ 
ture  which  contains  a  500  mb  short  wave  responsible  for  the 
rapid  surface  deepening  (Fetterssen  and  Smebye,  1971).  in 
the  type-A  case,  the  flow  aloft  is  predominantly  zonal 
through  the  early  explosive  deepening  phase  with  positive 
vorricity  advection  of  minor  significance.  Low  level  ther¬ 
mal  advection  plays  the  important  role  in  the  deepening 
process.  This  thesis  is  believed  to  be  one  of  the  first 
detailed  analyses  of  this  class  of  maritime  storms.  In 
addition,  the  vast  majority  of  explosive  deepeners  form  in 
the  North  Pacific  Ocean  (Sanders  and  Gyakum,  I960).  Inas¬ 
much  as  the  most  severe  weather  is  directly  associated  with 
these  severe  storms,  this  study  assumes  particular  signifi¬ 
cance  in  attempting  to  understand  this  meteorological 
phenomenon  in  order  to  mirimize  its  disastrous  effects. 

The  specific  approach  employed  in  this  thesis  is  the  use 
of  QLD  techniques  in  an  isobaric  spherical  coordinate  system 
which  is  translated  with  the  storm.  Horizontal  advection 
associated  with  wave  translation  is  isolated  from  the  rela¬ 
tive  advection  and  the  divergence  of  the  transport  associ¬ 
ated  with  cyclone  development.  The  fluctuations  of 
meteorological  variables  are  then  analyzed  through  the  use 
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of  a  cylindrical  budget  volume  centered  on  and  translating 
with  the  cyclone  surface  center.  Vertical  distributions, 
lateral  exchanges  and  sources  and  sinks  of  cyclone  proper¬ 
ties  resulting  from  purely  developmental  processes  are  then 
analyzed.  Application  of  QLD  techniques  in  studying  cyclone 
development  (Hash,  1978)  has  proven  valuable  in  diagnoses  cf 
poorly  forecast  storms,  and  Roman  (1981)  and  Conant  (1982) 
have  demonstrated  QLD  applicability  in  detailed  analyses  of 
east-coast  cyclogenesis. 

The  main  objectives  of  this  thesis  are: 

•  A  synoptic  investigation  of  the  classical  frontal  wave 
cyclogenesis  in  the  Ncrth  Pacific  Ocean; 

•  A  synoptic  assessment  of  the  validity  of  the  First  garp 
Global  Experiment  (FGGE)  data  and  the  quality  of  the 
FGGE  analyses  for  depicting  storm  dynamical  processes; 

•  A  diagnostic  investigation  of  storm  dynamics  through 
the  use  of  QLD  budget  analyses  of  mass  and  vorticity. 

The  detailed  analyses  of  these  processes  is  further 
aided  by  an  exceptionally  strong  supporting  data  base.  The 
case  under  study  occurred  during  the  first  special  observa¬ 
tion  period  (SOP-1)  of  the  First  GARP  Global  Experiment 
(FGGE).  This  FGGE  data  set  was  assimilated  by  the  European 
Center  for  Medium-range  Heather  Forecasting  (ECMWF)  and 
obtained  through  the  auspices  of  the  Goddard  Space  Flight 
Center.  This  data  set  is  special,  in  that  observations  from 
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aircraft,  ships,  buoys  and  satellites  gathered  during  an 
intense  collection  effort  were  Manually  checked  and  then 
systematically  integrated  into  a  dense  network  of 
observations.  Consequently,  the  supporting  data  base  ana¬ 
lyzed  in  this  thesis  effort  is  exceptionally  broad-based, 
complete,  and  integrated  into  a  well-diversified 
observational  network. 

Quantitative  analysis  cf  FGG E-type  data  bases  gives  spe¬ 
cific  insights  into  development  processes  inherent  in  mari¬ 
time  explosive  cyclogenesis.  In  the  following  chapters, 
examination  of  the  FGGE  data  (Chapter  2)  and  an  abbreviated 
description  cf  QLD  procedures  (Chapter  3)  are  presented. 
Chapter  4  discusses  the  cyclone  under  study  through  a  synop¬ 
tic  overview  which  highlights  significant  flow  patterns  and 
associated  upper  and  lower-level  features,  chapters  5  and  6 
present  the  results  of  the  budget  calculations  and  Chapter  7 
contains  specific  conclusions  from  these  budget  calculations 
and  discusses  recommendations  for  further  research. 
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II.  ENVIRONMENTAL  DAT  A  E  ASJ  DESCRIPTION 


&.  EGGE  DATA 

Curing  the  first  special  observation  period  (SOP-1)  of 
the  Global  Heather  Experiment,  the  largest  concentration  of 
meteorological  resources  ever  assembled  was  deployed  to 
observe  the  atmosphere  and  oceans.  The  quality  of  the  data 
generated  from  this  effort  marked  a  milestone  in  the 
development  of  numerical  weather  prediction  models  and 
synoptic  meteorology  in  general.  By  any  standard,  this 
experiment  shewed  an  unparalleled  example  of  scientific 
cooperation  and  excellent  resource  management. 

The  FGGE  Level  Ill-b  data  set  used  in  this  thesis  was 
constructed  from  observations  integrated  through  a  data 
assimilation  system  at  ECMWF.  Bengsston  et  al.,  (1982) 
describe  the  system  as  a  three-dimensional,  multivariate, 
optimum  interpolation  scheme  using  a  nonlinear  normal-mode 
initialization  and  an  automatic  system  for  data  checking. 
The  direct  field-level  observations  (Level  I-b)  were  trans¬ 
formed  into  the  basic  meteorological  variables  (Level  Il-b)  , 
which  constituted  the  basic  meteorological  product  of  the 
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experiment.  a  15- level  model  with  horizontal  resolution  of 
1.  875  degrees  of  latitude  and  longitude  was  used  to  trans¬ 
form  the  level  Il-b  data  to  Level  Ill-b.  At  the  Ill-b 
level,  the  observations  were  manually  checked  and  objec¬ 
tively  analyzed  into  a  set  of  dynamically  derived  meteorolo¬ 
gical  fields.  The  flowchart  of  this  data  assimilation 
process  is  depicted  in  Fig.  1.  The  FGGE  Level  Ill-b  data 
production  effort  was  completed  in  June  1981  and  the  data 
fields  were  catalogued  for  6-h  observations  for  all  standard 
levels  up  to  and  including  10  mb.  A  detailed  description  of 
this  assimilation  scheme  can  be  found  in  Bengsston  et  al., 
(1985)  . 

The  FGGE  data  base  is  currently  undergoing  validation  by 
a  variety  cf  groups.  Homan  (1981)  and  Conant  (1982)  demon¬ 
strated  its  validity  on  a  case  of  east-coast  cyclogenesis. 
In  an  independenx  assessment  of  the  FGGE  data  set,  Salem  et 
al.,  (1985)  present  a  detailed  analysis  in  connection  with 
the  coverage  and  accuracy  of  the  remote  sensing  portion  of 
the  network  covering  the  oceanic  regions.  Analyses  based  on 
these  inputs  provided  the  following  conclusions: 
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•  The  FGGE  observing  system  has  a  positive  effect  on 
forecast  skill  using  Si  criteria; 


•  Significant  corrections  to  the  300  mb  first  guess 
fields  cccur  over  oceanic  regions  with  the  FGGE  data; 


•  Nearly  all  of  the  forecasts  that  retain  useful  skill  by 
day  five  shew  a  positive  impact  from  the  FGGE  system; 


•  In  certain  synoptic  situations,  the  FGGE  observing  sys¬ 
tem  is  capable  of  measuring  the  atmosphere  Globally 
with  sufficient  accuracy  for  current  MNP  models  to  make 
useful  extended- range  forecasts  beyond  one  week. 


These  and  other  conclusions  as  well  as  a  detailed  descrip¬ 


tion  of  the  analysis  scheme  can  be  found  in  Halem  et  al.. 


(1981)  and  Baker  et  al. ,  (1981).  It  is  evident  from  these 


and  ether  studies  that  the  FGGE  Level  Ill-b  data  set  pro¬ 


vides  a  streeg  foundation  for  accomplishing  the  objectives 
of  this  thesis. 

A  summary  of  the  FGGE  data  availability  by  field-type 
and  level  is  depicted  in  Table  I.  The  advantages  of  this 


type  of  augmented  data  support  over  the  standard  data  base 


include; 


•  A  considerably  richer  data  base,  particularly  over 
data-sparse  oceanic  areas; 


•  The  high  resolution  grid  is  convenient  for  depicting 
the  smaller  synoptic-scale  features  associated  with 
cyclcgenesis; 


•  Basic  and  derived  fields  (vertical  velocity,  humidity, 
etc.)  are  included  in  the  data  set; 
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•  The  full  area  of  cyclogenesis  and  storm  translation 
enjoys  sinilar  data  coverage. 

These  significant  advantages  demonstrate  the  value  of  the 
FGGE  Level  Ill-b  data  set  in  quantitative  analyses  of  store 
features. 

B.  SUMMARY  OF  ANALYSES  AND  PROGNOSES  UTILIZED 

Three  sets  of  analyses  are  utilized  ir.  this  thesis: 

•  National  Meteorological  Center  (NMC)  analyses  and  prog¬ 
noses; 

•  Fleet  Numerical  Oceanography  center  (PNOC)  analyses  and 
prognoses; 

•  ECMWF  Level  m-fc  analyses. 

These  analyses  are  used  in  conjunction  with  the  synoptic 
discussion  (Chapter  4)  and  for  the  verification  of  the 
forecast  problem  alluded  tc  earlier.  A  complete  listing  of 
all  analyses  and  prognoses  utilized  in  this  study  is  con¬ 
tained  in  Table  n.  A  complete  description  of  the  FNOC 
products  can  be  found  in  the  U.s.  Naval  Heather  Service 
Numerical  Environmental  Products  Manual  (NAVAIR  50-1G-522). 
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III.  CDASI-LAGBANGIAN  DIAGNOSTIC  (QLDI 
TECHN1Q 


To  derive  generalized  budget  analysis  techniques  for 
processes  related  to  cyclcne  development,  changes  in  the 
budget  volume  properties  due  to  migration  of  the  system  must 
be  separated  tc  isolate  true  developmental  trends.  In  the 
past,  synoptic-scale  diagnostic  approaches  (a.g.  Rung  and 
Baker,  1975)  evaluated  storm  processes  relative  to  a  fixed 
coordinate  system.  This  type  of  approach  is  effective  for 
selected  time  periods  when  storm  translation  is  negligible, 
but  suffers  markedly  when  applied  to  rapidly  translating 
systems  of  appreciable  intensity.  In  particular,  changes 
may  cccur  due  tc  system  movement,  system  development,  or 
both.  The  only  accurate  approach  to  assess  true  developmen¬ 
tal  fluctuations  associated  with  cyclone  development  is  tc 
construct  a  moving  budget  volume  centered  on  some  distin¬ 
guishing  storm  feature.  In  this  way,  horizontal  advection 
associated  with  wave  translation  may  be  isolated.  In  addi¬ 
tion,  any  vertical  redistribution  of  a  budget  property  may 
be  analyzed  in  terms  of  vertical  transport  processes  within 
the  cyclone  budget  regime.  QLD  techniques  employed  in  this 
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thesis  provide  the  means  fcr  this  type  of  approach.  Reman 
(1981)  and  Conant  (1982)  have  demonstrat ad  the  validity  of 
this  type  cf  approach  on  a  casa  of  east-coast  cyclogenesis. 

To  study  the  fluctuations  in  any  storm  property,  a 
generalized  budget  equation  is  constructed  in  terms  of  the 
storm  budget  coordinate  system.  The  equation  consists  of 
the  time  rate  of  change  of  the  desired  budget  quantity 
(df/dt)  and  the  lateral  flax  transport,  vertical  flux 
transport  and  sources/sinks  (see  Table  III) .  The  equation 
is  applicable  to  the  storm  budget  volume  in  spherical  coor¬ 
dinates.  The  spherical  portion  of  the  atmosphere  incorpo¬ 
rated  within  the  budget  volume  is  approximated  by  a  cylinder 
fcr  ease  in  computations  and  definition  of  budget  proper¬ 
ties.  This  approximation  is  valid  for  the  relatively  thin 
slab  of  atmosphere  (1000-100  mb)  selected  for  budget 
analyses. 

The  formulation  of  the  storm  budget  coordinate  system 
and  details  cf  the  budget  equation  notation  are  discussed  by 
Wash  (1978).  The  storm  budget  volume  is  defined  in  terms  of 
10  mandatory  pressure  levels  in  the  vertical,  and  in  the 
horizontal  by  outwardly  directed  radii  originating  from  a 
central  point.  For  example,  the  budget  property  at  a  point 
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may  te  analyzed  ty  describing  the  level/layer  (1000/925  mb) , 
the  budget  radius  in  degrees  of  latitude  from  the  storm 
center  (6)  ,  and  the  quadrant  desired  (grid  points  1-36) . 
The  mass  and  vorticity  budget  equations  employed  in  this 
analysis  are  depicted  in  Tables  IV  and  V  respectively. 

In  the  mass  budget  (Table  IV) ,  the  general  equation  rep¬ 
resenting  the  change  in  mass  with  time  is  composed  of  con¬ 
tributions  from  lateral  and  vertical  transports.  The 
lateral  transport  term  represents  the  net  horizontal  conver¬ 
gence  of  mass  flux  within  the  budget  volume.  The  vertical 
transport  is  calculated  through  the  use  of  derived  vertical 
velocity  fields  and  is  recessarily  a  reflection  of  net 
lateral  flux  into  the  volume.  Application  cf  boundary  con¬ 
ditions  that  vertical  motion  is  zero  at  the  top  and  bottom 
surfaces  of  the  budget  volume  simplifies  the  equation,  as 
shewn  in  the  middle  of  Table  IV. 

Several  useful  relationships  can  be  derived  from  this 
simplified  equation  and  the  boundary  conditions.  Contribu¬ 
tions  to  the  mass  time  tendency  in  iso  baric  coordinates 
arises  chiefly  from  pressure  fluctuations  in  the  surface 
layer,  as  the  ether  layers  in  the  budget  are  fixed  pressure 
surfaces.  This  term  is  computed  easily  and  with 
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considerable  accuracy.  The  net  lateral  mass  transports  in 
the  layers  are  calculated  from  the  wind  fields  and  necessar¬ 
ily  balance  the  vertical  transports  through  continuity. 
Errors  in  the  wind  fields  cause  the  vertically  integrated 
lateral  mass  transport  to  differ  from  the  vertically 
integrated  mass  tendency,  and  cause  a  residual  in  the 
budget. 

Application  of  the  boundary  conditions  to  the  divergence 
of  the  vertical  transport  causes  the  vertical  sum  to  vanish, 
and  thereby  isolates  the  residual  correction  which  must  be 
applied  layer  by  layer  tc  force  the  vertical  sum  of  the 
divergence  tc  zero.  The  O'Brien  (1970)  linear  correction 
scheme  is  utilized  to  correct  for  these  residual  imbalances. 
The  upper  layers  are  weighted  more  due  to  the  greater  inac¬ 
curacies  of  the  wind  measurements  at  these  levels.  These 
corrections  are  applied  tc  the  layer  mass  transports,  and 
are  combined  with  the  translation  correction  (discussed 
earlier)  to  give  the  final  corrected  form  of  the  mass 
transport  term.  The  vertical  velocities  representing  the 
vertical  mass  transport  are  then  recovered  through 
continuity. 
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In  the  vorticity  budget  formulation,  the  contributions 
to  vorticity  changes  occur  from  lateral  and  vertical  fluxes, 
sources  and  sicks  (Table  V,  Section  A)  .  The  flux 
(transport)  terms  can  be  expanded  into  divergence  and  advec- 
tion  contributions  by  the  use  of  vector  identities  (Table  7, 
Section  B)  ,  which  facilitates  further  insight  into  the  phys¬ 
ical  processes  they  represent.  The  source/sink  terms  arise 
from  net  generation/dissipation  of  vorticity  within  the 
budget  volume,  and  in  iscbaric  coordinates  are  comprised  of 
divergence,  tilting  and  friction  components.  Notice  that 
three  divergence  terms  arise  from  this  expansion  of  the 
general  vorticity  budget  eguation.  The  horizontal  transport 
component  is  identical  to  the  source  divergence  term,  and 
the  horizcntal  and  vertical  transport  components  are  equal 
and  opposite.  Thus,  the  net  contribution  of  divergence  in 
the  budget  is  manifest  in  the  source  term.  A  specific 
breakdown  of  the  partitioned  vorticity  budget  equation  is 
presented  in  Section  B  of  Table  V. 

To  isolate  the  contributions  of  both  the  rotational  and 
irrotational  wind  components  to  the  transport  of  a  desired 
budget  quantity,  mean  and  eddy  modes  (analogous  to 
divergence  and  advection  terms  ir.  the  vorticity  equation)  , 
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are  derived  from  the  transport  terms  (Table  VI)  .  The  pur¬ 
pose  of  this  technique  is  to  isolate  the  changes  in  the  mean 
budget  properties  from  these  occurring  on  the  eddy  scale. 
Thus,  a  comprehensive  analysis  of  the  lateral  and  vertical 
flux  terms  can  be  made  by  separating  them  into  mean  and  eddy 
mode  contributions  and  comparing  them  with  the  divergence 
and  advection  terms  for  analysis  of  the  physical  processes 
they  represent.  The  relation  of  the  mean/eddy  mode  and 
divergence/aavection  partitions  to  the  total  transport  is 
also  illustrated  in  Section  B  of  Table  V. 

The  budget  analysis  begins  with  the  choice  of  the  par¬ 
ticular  property  to  be  analyzed.  The  next  step  in  the 
budget  calculations  consists  of  interpolating  the  u  and  v 
components  of  the  PGGE  data  set  onto  the  cylindrical  budget 
volume  grid.  Next,  the  interpolated  winds  are  converted  to 
tangential  and  radial  components  relative  to  the  outer  vol¬ 
ume  boundary  and  corrected  for  storm  motion.  In  this  form, 
the  data  are  ready  to  be  used  in  the  budget  equations 
discussed  above.  Derived  fields  (moisture,  vertical  veloc¬ 
ity,  kinetic  energy,  etc.),  not  directly  used  in  budget  com¬ 
putations  can  be  used  tc  supplement  insight  into  storm 
dynamics.  Interesting  comparisons  between  FGGE  derived 
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vertical  velocity  and  QLD  calculated  vertical  velocity,  for 


example,  serve  to  validate  the  QLD  technique  approach. 


The  QLD  approach  provides  insight  into  the  following  key 


properties : 


•  The  Magnitude  and  significance  of  the  horizontal  mass 
transport  structure; 


•  The  magnitude  and  significance  of  the  horizontal  vor- 
ticity  transport  structure; 


•  The  vertical  transport  structure  of  mass  and  vorticity 
derived  from  vertical  motion  estimates; 


•  The  magnitude  and  nature  of  source/sink  contributions 
to  budget  analyses. 


The  next  chapter  discusses  the  development  sequence  of  the 


storm  thrcugh  a  synoptic  discussion  utilizing  synoptic 


fields,  analyses  and  forecasts,  and  satellite  imagery. 


28 


IV.  SYNOPTIC  OVERVIEW 

A.  GENERAI 

Oceanic  cyclcgenesis  mecbani sms  can  in  part  ba  portrayed 
by  graphical  representation  of  pertinent  meteorological  and 
oceanic  fields.  This  chapter  explores  the  synoptic  condi¬ 
tions  associated  with  a  specific  case  of  North  Pacific 
explosive  cyclcgenesis  during  12-17  January  1979.  An 
attempt  is  made  to  gain  insight  into  formative  processes  and 
possible  factors  related  to  the  serious  shortcomings  in  the 
model  repesentaticn  of  the  dynamics  associated  with  these 
vigorous  events.  Use  is  made  of  FGGE  data  analyses.  Defense 
Military  Satellite  Program  (DMSP)  imagery  and  FNOC  sea 
surface  temperature  (SST)  analyses.  Comparisons  are  made 
between  FGGE  ,  FNOC  and  NMC  sea  level  pressure  analyses  and 
prognoses  to  highlight  differences  in  forecast  deepening 
rates  and  analyzed  and  forecast  storm  positions.  Explosive 
deepening  criteria  were  adopted  from  Sanders  and  Gyakum 
(1980),  who  addressed  the  NMC  7-layer  PE  model  performance 
associated  with  "explosive"  prediction.  Predicted  deepening 
rates  for  the  cyclones  they  studied  were  consistently  less 


than  1/3  of  the  observed  rates,  and  North  Pacific  explosive 
deepening  east  of  the  dateline  was  virtually  missed  by  the 
model.  An  assessment  of  the  FNOC  5-layer  PE  model  perform¬ 
ance  (Syakua,  1981)  demonstrated  similar  results. 
Specifically,  the  model  captured  only  1/3  of  the  observed 
central  pressure  falls  for  all  time  periods.  In  addition, 
the  most  accurate  of  all  the  forecasts  were  shown  to  be 
valid  between  24-36  h.  In  comparing  the  FGGE  analyses  with 
the  FNOC  and  NMC  analyses,  it  should  be  noted  that  the  FNOC 
analyses  and  prognoses  are  derived  from  the  5-layer  PE 
coarse  mesh  (381  tan  hcrizcntal  resolution)  model.  In  con¬ 
trast,  the  NHC  analyses  are  strictly  hand-analyzed  Northern 
Hemisphere  strips,  while  the  prognoses  are  derived  from  the 
7-layer  PE  fine  mesh  (190.5  km  horizontal  resolution)  model. 

3.  MODEL  PREDICTION  PERFCRMANCS 

The  storm  track  (at  12-h  increments)  is  presented  in 
Fig.  2a.  The  cyclone  develops  southeast  of  Japan  in  the 
vicinity  30N,  135E  and  transits  toward  a  final  position  off 
the  west  coast  of  Alaska.  Notice  the  significant  transla¬ 
tion  increase  after  0000  GMT  13  January.  The  5-day  evolu¬ 
tion  (0000  GMT  12  January  1979  to  0000  GMT  17  January  1979) 
in  central  sea-level  pressure,  as  analyzed  from  the  FGGE 
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2b. 


data,  SMC  and  FNOC  analyses,  is  illustrated  in  Fig. 
Assuming  the  representativeness  of  the  FGGE  data,  both 
series  of  FNOC  and  NMC  analyses  fail  markedly  in  represent¬ 
ing  the  intense  deepening  phase  (0000  GMT  13  to  0000  GMT  15 
January  1919)  .  There  is  close  agreement  between  FGGE  and 
NHC  analyses  at  other  times.  It  appears  that  specific 
inputs  to  the  FGGE  data  set  (mentioned  earlier)  not 
available  to  FNOC  or  NMC  for  real-time  analyses  provide  a 
more  accurate  portrayal  of  the  deepening  period.  Close 
inspection  of  the  observations  included  in  the  FGGE  analyses 
shows  that  the  deepening  phase  is  well  documented  at  all 
levels.  At  1200  GMT  1 «  January  1979,  for  example,  as  many 
as  seven  sea-level  observations  from  land  and  sea  stations 
ar 3  available  in  the  north  and  east  quadrants  of  the  stcrm 
and  form  the  surface  data  network.  In  the  upper  levels, 
satellite-derived  temperature  and  height  anomalies,  and  250 
mb  aircraft  winds  and  temperatures  combine  to  form  a  dense 
observational  network.  In  the  NMC  analysis,  ship  reports  in 
the  area  were  conspicuously  absent,  as  the  analysis  appeared 
to  be  based  upon  land  stations  only.  The  FNOC  analysis, 
however,  contained  four  of  the  seven  FGGE  surface 
observations,  but  did  not  contain  the  two  observations 
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closest  to  tfce  stcrm  center,  and  this  caused  a  misrepresen¬ 
tation  of  the  surface  central  pressure.  Thus,  the  mors 
accurate  representation  of  the  deepening  phase  in  the  FGGE 
analysis  is  directly  attributed  to  a  superior  data  base  and 
not  just  tc  the  benefit  of  a  "better"  first  guess. 

A  comparison  between  the  36-h  sea-level  pressure 
forecasts  is  shown  in  Fig.  3a.  All  curves  show  the  general 
deepening  trends,  but  the  FNOC  and  NMC  forecasts 
consistently  lag  the  FGGE  analyzed  pressure,  and  are  highly 
variable  in  predicting  the  magnitude  of  the  central  pressure 
falls.  From  these  displays,  it  is  evident  that  both  models 
fail  significantly  in  the  timely  and  accurate  prediction  of 
stori  intensification,  especially  during  the  explosive 
stage.  This  example  illustrates  further  the  forecast  prob¬ 
lems  discussed  by  Sanders  and  Gyakum  (1980)  alluded  to 
ea  rlier. 

The  location  error  in  36-h  storm  positions  for  NMC  and 
FNOC  predictions,  relative  to  the  FGGE  analyzed  position,  is 
presented  in  Fig.  3b.  It  is  clearly  seen  that  significant 
errors  exist  in  both  forecasts  at  various  times  during  the 
deepening  phase.  Such  significant  errors  in  the  analyzed 
and  forecast  features  of  this  system  reflect  an  inability  of 
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current  operational  aodels  to  accurately  forecast  these 
events. 

C.  SYNOPTIC  DISCUSSION 

The  synoptic  analysis  approach  in  this  section  consists 
of  an  examination  of  the  basic  meteorological  fields  in 
three  displays  of  two  fields  each.  Display  A  contains  250 
mb  winds  and  isotachs,  B  presents  500  mb  heights  and  abso¬ 
lute  vorticity,  and  C  shows  1000-500  mb  thickness  and  sea- 
level  pressure  patterns.  This  approach  generally  allows 
only  subjective  insighns  into  features  of  interest  during 
the  storm  life  cycle.  The  quantitative  analysis  of  the 
dynamical  processes  is  presented  in  the  following  chapter. 
The  intent  of  this  section  is  to  trace  the  development  of 
the  system  through  a  synoptic  discussion  of  major  features 
contriouting  to  overall  storm  development. 

Examination  of  atmospheric  conditions  at  0000  G?1T  12 
January  1979  in  the  area  southeast  of  Japan  shows  marked 
cyclcgenetic  potential.  The  surface  picture  (Pig.  4c)  is 
dominated  ty  a  broad  trough  (vicinity  30M,  14QE)  extending 
northeast-southwest  through  the  center  of  a  checkerboard 
pattern  of  fairly  strong  low  and  high  centers.  The  signifi¬ 
cant  features  include  the  Siberian  Hign  near  55N,  100E  (jus t 
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off  the  western  edge  of  the  chart)  which  is  very  strong  and 
stationary.  The  parent  lew  in  the  vicinity  of  53N,  17QE  is 
now  guasi-stationary  and  filling  (950  ah).  A  large  maritime 
high  (1030  mb)  dominates  the  north  central  Pacific  (vicinity 
30  N,  178W)  ,  and  the  northern  fringe  of  Typhoon  Alice  is  just 
visible  near  20N,  140E)  .  The  combination  of  flow  from  these 
systems  serves  to  funnel  cold,  polar  continental  air  south¬ 
ward  towards  northward  flcwing  warm,  moist  tropical  air, 
which  enhances  the  frontogenes is  potential.  A  general 
destabilization  of  the  polar  air  occurs  north  of  the  area  of 
strong  low  level  baroclinicity  shown  by  the  concentration  of 
thickness  contours  just  east  of  Japan. 

In  the  mid-troposphere  (Fig.  4b)  ,  two  distinctive  fea¬ 
tures  are  evident.  The  flow  over  the  western  Pacific  is 
dominated  by  an  extensive  region  cf  low  heights  associated 
with  a  guasi-stationary  long  wave  trough,  which  is  a  portion 
of  a  four-wave  planetary  pattern.  This  significant  feature 
accounts  for  the  steady  progression  of  mobile  short  waves 
and  associated  surface  features  through  the  region  during 
the  preceding  week,  in  addition  to  supporting  the  steady 
flow  of  cold  air  off  the  Siberian  continent.  The  second 
major  feature  of  interest  is  the  predominantly  zonal  flew 
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over  tha  entire  incipient  region  with  areas  of  weak  and 
small  absolute  vcrticity  upstream. 

In  the  upper  troposphere  (Fig.  4a),  the  presence  of  a 
broad  jet  (  >  50  m/sec)  aloft  is  a  further  manifestation  of 
the  fcaroclinicity  of  the  lower  troposphere.  As  yet,  signif¬ 
icant  jet  streak  (  >  75  m/sec)  propagation  into  the  region 
remains  well  upstream,  and  therefore  is  not  a  factor  in  the 
low-level  development. 

Examination  of  the  visual  (Fig.  5)  and  infrared  (Fig.  6) 
satellite  imagery  valid  at  0041  GBT  12  January  1979  shows  a 
couplet  of  similar  low  level  cloud  masses  (labelled  A  and  B) 
southeast  cf  Japan.  Lack  of  significant  vertical 
development  is  evident  on  the  infrared  imagery,  in  contrast 
to  the  area  cf  enhanced  convection  (labelled  C)  associated 
with  an  area  of  positive  vcrticity  advection  (PVA)  upstream 
over  Korea  (Fig.  4b) .  Cloud  mass  B  is  the  precursor  cf  a 
weak  cyclone  which  moves  to  the  east  and  rapidly  dissipates. 
Cloud  mass  A  is  the  precursor  of  the  cyclone  analyzed  in 
this  study.  The  evolving  low-amplitude  "s-shape"  configura¬ 
tion  with  a  distinct  edge  on  the  western  flank  is  indicative 
of  increasing  thermal  and  low  level  absolute  vorticity  gra¬ 
dients  with  time  (Weldon,  1977).  Although  fairly  small,  the 
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cloud  mass  is  beginning  tc  tales  on  the  familiar  "leaf-shape" 
which  Weldon  (1977)  characterizes  as  typical  of  incipient 
oceanic  cyclogenesis. 

During  the  nert  24  h  ending  at  0000  GHT  13  January  1979, 
surface  development  continues  (vicinity  35N,  15GE)  into  a 
closed  center  (995  mb).  The  incipient  system  (Fig.  7c) 
appears  as  a  broad  area  of  low  pressure,  and  the  development 
of  a  weak  thermal  ridge  near  the  surface  center  indicates 
frontogenesis  continues  during  the  period.  Ships  in  the 
vicinity  reported  seas  on  the  order  of  15  ft,  with  periods 
ranging  between  3-5  s,  which  indicates  a  duration  limited 
wave  spectrum.  Of  interest  is  the  continued  strengthening 
of  the  Siberian  high  (now  1080  mb)  which  steadily  pumps  cold 
air  southward,  in  association  with  the  backside  circulation 
of  the  stationary  c^'.d  low  to  the  east.  Destabilization  of 
the  surface  layer  continues  as  cold  air  streams  southward 
over  the  warm  ocean  surface,  and  a  polar  low  has  become 
organized  ever  the  southeast  tip  of  Japan  due  west  of  the 
incipient  system  (indicated  on  detailed  NMC  analyses,  not 
shown  hare) . 

In  the  mid-troposphere  (Fig.  7b|  ,  the  500  mb  low 
anchored  off  the  Kamchatka  Peninsula  has  deepened.  A  short 
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wave  with  at  associated  broad  absolute  vorticity  area  has 
formed  over  Japan,  approximately  800  nm  upstream  of  the 
incipient  system.  A  smaller,  less  distinct  vorticity  area 
has  formed  directly  over  the  surface  system. 

Increased  vorticity  in  the  mid-troposphere  can  also  be 
linked  to  the  approach  of  a  jet  streak  from  the  west.  At 
this  stage,  however,  the  position  of  the  jet  streak  indi¬ 
cates  only  minimal  influence  on  the  low-level  circulation, 
as  the  jet  core  is  directly  superposed  in  the  vertical,  and 
no  clear  pattern  of  cyclonic  shear  or  diffluence  can  be  seen 
at  250  mb  (Fig.  7a)  . 

The  ircrsased  vertical  development  over  the  period  can 
be  seen  in  the  DMSP  satellite  imagery  (Figs.  8  and  9)  valid 
at  2341  GUT  12  January  1979.  The  developing  system  is 
labeled  A  and  is  centered  approximately  600  nm  southeast  of 
Japan.  Cold  plumes  extend  for  virtually  the  entire  breadth 
in  the  infrared  image  (Fig.  9),  with  cloud  tops  becoming 
progressively  warmer  toward  the  tapered  end  to  the  south¬ 
west.  Enhanced  convection  and  continued  development  of  the 
leaf  structure  is  shown  in  the  visual  image  (Fig.  8).  The 
region  under  the  leaf  appears  to  contain  marked  oarcclinic- 
it  y  in  a  zone  in  advance  of  the  propagating  ■*<?-  streak 
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aloft,  which  is  similar  tc  the  conditions  observed  by  Weldon 
(1  977)  in  a  case  study. 

Curing  the  following  24-h  period  ending  0000  GMT  14  Jan¬ 
uary  1979,  continued  development  is  evident  in  the  low  and 
middle  troposphere.  Central  surface  pressure  (Fig.  10c) 
falls  steadily  tc  977  mb,  which  is  a  decrease  of  19  mb.  If 
one  adopts  the  explosive  criteria  of  Sanders  and  Gyakum 
(1980)  for  a  latitude  of  40  degrees,  the  storm  has  entered 
the  explosive  stage  during  this  period.  The  surface  low 
moves  rapidly  toward  the  northeast,  to  a  position  close  to 
the  maximum  SSI  gradient  (Fig.  11).  Close  proximity  to  this 
gradient  fulfills  another  observation  of  explosive 
development  described  ty  Sanders  and  Gyakum  (1980). 
Analysis  of  the  thickress  contours  (Fig.  10c)  indicates  a 
pattern  of  diffluence  10  degrees  of  longitude  downstream 
from  the  surface  center.  The  diffluence  at  10  degrees 
upstream  is  ccnsiderably  greater.  Bjerknes  (1954)  shewed 
that  these  conditions  favor  formation  of  a  diffluent  trough 
in  the  upper  levels  within  the  cyclonic.-  shear  zone  north  of 
the  jet  streak  axis.  The  NW-SE  tilt  of  the  500  mb  trough 
evident  by  1200  GHT  14  January  (Fig.  16b) ,  combined  with  the 
favorable  superposition  of  the  jet  streak  support  the 
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diffluent  trough  formation  argument.  In  addition,  a  combi¬ 
nation  of  the  superposition  of  the  upper  level  divergence 
zone  ever  the  surface  center  and  the  approach  of  the  upper- 
level  thermal  trough  in  close  proximity  to  the  500  mb  trough 
promotes  self-amplification  through  the  mid-troposphere. 
Mid-level  vorticity  increases  due  to  the  movement  of  concen¬ 
trated  cyclonic  shear  vorticity  into  the  budget  volume  as 
the  jet  enters  the  diffluent  trough.  As  a  result,  increased 
vertical  development  occurred  from  vorticity  increases  as  a 
result  of  jet  streak  interaction  and  self-amplif ication . 

Indeed,  the  first  evidence  of  mid-tropospheric 
development  is  visible  in  Fig.  10b.  The  vorticity  pattern 
has  an  expanded  area  of  significant  vorticity  associated 
with  the  mid-tropospheric  circulation.  A  region  of  strong 
vorticity  also  exists  upstream  northeast  of  Japan,  but  that 
center  appears  tc  be  associated  with  the  polar  low  forming 
to  the  west  (Fig.  10c),  as  indicated  earlier.  The  interest¬ 
ing  feature  tc  notice  is  the  absence  of  the  short  wav®  that 
was  present  ever  Japan  24  h  earlier,  which  would  seem  to 
indicate  that  further  explosive  development  from  the  short 
wave  PVA  alcft  is  unlikely. 
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At  250  mb  (Fig.  10a),  however,  the  arrival  of  the  jet 
streak  (maximum  velocity  85  m/sec)  in  favorable  superposi¬ 
tion  over  the  surface  system  acts  to  promote  the  lower-level 
circulation.  Diffluence  in  the  region  of  the  jet  stream 
split  and  cyclonic  wind  shear  to  the  left  of  the  jet  axis 
combine  to  induce  positive  circulation  increases  and  strong 
vertical  motion  through  the  mid-troposphere.  Increased 
vertical  motion  enhances  vertical  redistributions  of  all 
properties,  including  an  increase  in  the  convergence  in  the 
surface  layer,  and  hence,  explosive  development  continues. 

Inspection  of  satellite  imagery  (Figs.  12  and  13  valid 
1405  GMT  13  January  1979  and  Figs.  14  and  15  valid  at  0023 
GMT  14  January  1979)  shows  a  discernable  comma-shape 
formation.  In  Figs.  12  and  13,  the  southern  boundary  of  the 
jet  axis  is  delineated  by  subsidence  along  the  cirrus  band 
paralleling  the  lew- level  front.  The  jet  crosses  just  south 
of  the  vortex  at  A.  The  very  ccld  tops  evident  near  the 
vortex  center  and  along  the  frontal  band  indicate  strong 
convection  throughout  the  region.  An  interesting  feature  in 
Figs.  14  and  15  is  the  relatively  low  level  clouds  (labeled 
A)  that  appear  as  fingers  extending  from  the  vortex  center 
and  look  similar  to  feeder  bands.  When  coupled  with  the 
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extensive  cirrus  shield  aloft,  the  cloud  pattern  presents  a 
hurricane-like  appearance.  Two  significant  changes  have 
occurred  during  the  past  10  h.  The  area  of  deep  convection 
has  expanded  markedly  both  near  the  vortex  and  along  the 
frontal  band,  which  indicates  continued  strengthening. 
Second,  the  storm  has  organized  to  the  point  that  clear 
slots  have  formed  in  the  cold  air  behind  the  front,  which 
indicates  the  penetration  of  cold  and  subsiding  air  toward 
the  vertex. 

Curing  the  next  12-h  period  ending  1200  GMT  14  January 
1979,  the  surface  pressure  falls  an  additional  11  mb  to  965 
mb  (Fig.  16c).  Notice  the  shift  of  the  reference  grid  10 
degrees  to  the  nerth  and  30  degrees  to  the  east  in  this  and 
subsequent  figures.  The  increased  amplitude  of  the 
thickness  ridge  implies  increased  circulation  in  the  lower 
troposphere,  most  probably  as  a  result  of  warm  advection 
and  extreme  latent  heating  from  de«p  convection.  Sensible 
heat  input  presumably  weakens  during  the  period  as  the  storm 
center  moves  north  of  the  maximum  5ST  gradient  (shewn 
earlier  in  Fig.  1 1)  . 

At  500  mb  (Fig.  16  b)  ,  self-amplification  continues  with 
the  development  of  the  upper-level  trough  and  associated 
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vorticity  maximum.  The  apparent  lack  of  appreciable  vortic- 
ity  advection  at  this  level  suggests  that  continued  intensi¬ 
fication  of  the  system  must  be  manifest  in  either  the  lower 
or  upper  levels.  The  most  likely  factors  which  serve  to 
maintain  the  deepening  phase  consist  of  a  combination  of 
processes:  jet  stream  interaction  at  250  mb;  and  the 
extreme  surface  convergence  generating  pockets  of  sustained 
deep  convection. 

An  increase  in  cyclonic  vorticity  appears  at  the  250  mb 
level  during  this  time  frame  (Pig.  16a).  Cyclonic  curvature 
and  shear  to  the  left  of  the  jet  axis  continues  to  promote 
the  circulation  discussed  above.  The  wind  speeds  in  the  jet 
have  diminished  slightly  in  intensity  due  to  the  decreased 
baroclini city  through  the  lower  troposphere. 

Strong  development  of  the  system  during  the  previous  12 
h  is  indicated  in  Figs.  17  and  18  (valid  at  1205  GST  14  Jan¬ 
uary  1979)  .  Vertical  velocities  within  the  storm  volume 
increase  during  the  period,  as  deep  convection  becomes 
intense  and  widespread  surrounding  the  vortex  and  along  the 
entire  frontal  band.  The  growth  of  convective  activity 
indicates  that  development  continues  upward  through  the 
troposphere  during  the  period,  as  the  jet  crosses  the 
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frontal  band  at  a  more  northerly  point.  In  addition,  the 
"hook"  or  tip  of  the  cloud  band  in  the  northwest  sector 
develops  markedly.  In  typical  cases  (Weldon,  1977)  ,  these 
conditions  indicate  the  onset  of  the  nature  stage. 

Interestingly,  the  next  6-h  period  ending  1800  GUT  14 
January  1979  narks  the  most  explosive  period  cf  the  storm's 
development  cycle.  The  surface  pressure  (Fig.  19c)  falls  18 
mb  (averaged  between  FGGE  and  NMC  analyses).  Ship 
observations  at  1200  and  1800  GMT  were  limited  in  both 
analyses.  It  is  possible  that  seme  of  the  extreme  deepening 
observed  during  this  6-h  period  actually  occurred  before 
1200  GMT,  but  there  were  inadequate  ship  observations  to 
detect  it.  Thus,  the  magnitude  of  the  surface  pressure  fall 
may  be  less  than  mentioned  above  or  presented  in  Fig.  3b. 
In  ar.y  event,  the  nature  of  the  explosive  deepening  has 
caused  the  thickness  ridge  to  double  in  amplitude  along  an 
east-west  orientation,  indicating  that  cold  air  has  begun  to 
wrap  around  the  vortex  at  the  surface.  In  the  mid- 
tropesphere  (Fig.  19b),  the  positive  vorticity  center  elon¬ 
gates  along  the  flew  and  intensifies  slightly.  The 
elongation  of  the  vorticity  center  and  the  continued  digging 
of  the  mid-level  trough  indicate  the  system  should  weaken  in 
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■the  near  future  (Weldon,  1977).  At  250  mb  (Pig.  19a)  ,  evi¬ 
dence  still  suggests  that  the  jet  stream  is  in  favorable 
superposition  over  the  surface  system  and  is  continuing  to 
induce  vertical  lotion  and  associated  low-level  convergence. 
Although  the  jet  remains  fairly  stationary  and  strong,  the 
surface  lew  continues  northeast  at  a  high  speed,  which 
suggests  decreased  jet  influence  in  future  development. 

A  look  at  the  satellite  imagery  (Pigs.  20  and  21)  valid 
at  2124  GHT  14  January  1979  shows  the  storm  at  maximum 
development.  The  dry  tengue  aloft  has  penetrated  to  the 
vortex  center.  A  solid  deck  of  st ratocumulus  to  the  south¬ 
east  of  the  vortex  center  exists  in  the  strong  cyclonic 
surface  flow.  An  interesting  "hole"  appears  to  the  east  cf 
the  vertex  in  the  dense  frontal  band  clouds,  which  is  not 
due  to  orographic  effects  as  the  area  is  still  south  of  the 
Aleutian  chain.  A  detailed  look  at  the  preceding  satellite 
picture  (Fig.  17a)  shows  a  clear  slot  formed  in  about  the 
same  relative  pcsiticn.  Thus,  the  feature  is  associated 
with  the  rotation  of  the  cloud  band  around  this  dry  tongue 
in  the  streng  cyclonic  flew  at  the  surface,  and  gives  the 
appearance  of  the  formation  of  a  secondary  vortex  in  the 
frontal  band.  Further  inspection  of  this  imagery  shows  the 
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continuance  cf  strong  vertical  notion  induced  by  deep  con¬ 
vection  alcng  the  frontal  band,  and  to  the  north  and  east  of 
the  vortex  center.  The  penetration  of  cold  air  into  the 
vortex  is  visible  and  signals  the  start  of  the  occlusion 
phase. 

At  0000  GdT  15  January  1979,  the  surface  circulation  is 
dominated  by  the  surface  lew  over  the  Aleutian  Islands  (Fig. 
22c)  with  a  central  pressure  at  947  mb.  The  explosive 
period  has  ended  with  the  onset  of  the  occlusion  process. 
The  thickness  gradient  has  diminished  markedly,  and  a  polar 
low  is  organized  in  the  cold  air  to  the  west,  while  a  new 
disturbance  forms  alcng  the  baroclinic  zone  to  the  south¬ 
west.  At  500  mb  (Fig.  22b) ,  a  closed  circulation  has 
formed,  which  indicates  cold  air  influx  through  most  of  the 
troposphere.  A  broad  positive  vcrticity  pattern  exists,  but 
negative  vorticity  advection  into  the  region  is  inducing 
subsidence  and  filling.  The  system  is  vertically  stacked, 
which  also  indicates  the  fully  mature  stage.  At  250  mb 
(Fig.  22a),  development  cf  the  trough  is  evident  from  the 
cyclonic  rotation  of  the  wind  vectors.  The  jet  has  dimin¬ 
ished  in  intensity  due  to  the  reduced  low-level  baroclinic- 
ity  and  the  migration  of  the  surface  system  north  of  the  jet 
axis . 
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During  the  next  24  hcurs  ending  0000  GMT  16  January 
1979,  the  occlusion  process  continues.  The  surface  system 
(Fig.  23c)  has  stalled  over  the  west  Alaskan  coast  and 
filled  to  955  oh.  Ship  reports  in  the  vicinity  show  35  ft 
seas  with  30  ft  swell  and  indicate  that  the  storm  is  still  a 
potent  force.  Other  signs  indicate,  however,  that  the  storm 
is  beginning  to  dissipate.  The  decrease  in  the  thickness 
gradient  indicates  that  the  lower  tropospheric  front  is 
occluded.  The  500  mb  absclute  vorticity  (Pig.  23b)  decrease 
indicates  a  large  reduction  in  relative  vorticity  from  neg¬ 
ative  vorticity  advection  and  the  discontinuance  of  the  jet 
stream  influence  (Fig.  23a).  Extreme  height  falls  are 
reflected  in  the  250  mb  flow  as  the  wind  field  becomes 
closed,  indicating  an  upward  cyclonic  development  through 
this  level. 

Inspection  of  the  satellite  imagery  for  this  period 
(Figs.  24  and  25)  valid  at  1005  GMT  15  January  1979  shows  a 
fully  mature  occluded  system.  The  penetration  of  cold  air 
to  the  vortex  center  is  evident  along  with  a  marked  decrease 
in  deep  convective  activity  near  the  vortex  and  along  the 
frontal  band,  except  in  the  northern  fringe  where  the  strong 
cyclonic  flow  is  still  producing  some  low-level  convergence 
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and  associated  deep  convection.  The  developing  low  to  the 
south  is  just  visible  in  the  imagery. 

By  0000  GMT  17  January  1979,  the  surface  system  has 
filled  to  962  mb  and  the  thickness  gradient  is  small  (Fig. 
26c) .  At  500  mb  (Pig.  26b)  ,  the  low  has  combined  with  a 
stationary  feature  to  the  west  to  produce  a  vast  area  of  low 
heights.  Positive  vorticity  has  decreased  as  the  vertical 
gradients  have  deteriorated  markedly  while  the  system  slowly 
dies.  At  250  mb  (Fig.  26a),  the  increased  closed 
circulation  indicates  vertical  stacking  through  the  whole 
troposphere.  These  features  suggest  that  the  storm  has 
reached  the  dissipating  phase,  which  marks  the  termination 
of  the  syncptic  portion  of  the  analysis. 

It  is  interesting  to  notice  that  during  the  period  of 
lowest  central  pressure  (0000  GMT  15  January  1979)  ,  ship 
reports  in  the  region  indicate  40  ft  waves,  20  ft  swell  and 
38  kt  winds.  The  dramatic  deepening  of  this  and  similar 
storms  is  net  without  its  consequences.  These  storms  pose 
significant  threats  to  ships  and  land  facilities  within 
their  paths.  The  intent  of  the  next  chapter  is  to  pursue 
those  processes  qualitatively  described  in  this  chapter  by 
mass  and  vcrticity  budget  analyses.  The  budget  framework 
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provides  a  ccaputational  base  through  which  significant 
aspects  of  the  storm  are  explored. 


7 


MASS  BODGE?  ANALYSIS 


A.  GENERAL 

It  is  clear  from  the  synoptic  discussion  (Chapter  4) 
that  the  key  factors  necessary  for  cyclone  development  are 
inherently  related  to  dynamical  processes  occurring  in  the 
lower  tropospheric  regime.  One  such  vital  process  is  the 
destafcilizaticn  cf  the  surface  layer.  Another  includes  the 
horizontal  and  vertical  mass  circulation  structure,  as  pres¬ 
sure  falls  and  resultant  storm  intensity  are  directly  depen¬ 
dent  upon  the  nature  and  magnitude  of  the  mass  flux  into  and 
out  cf  the  budget  volume. 

The  mass  budget  portion  of  this  chapter  includes  a  quan¬ 
titative  approach  to  determine  the  destabilization  and  hori¬ 
zontal  and  vertical  mass  transport  within  the  budget  volume. 
The  stability  analysis  is  based  on  area-averaged  potential 
temperatures  fcr  the  layer.  For  example,  the  layer  labelled 
775  mb  represents  the  850-700  mb  layer.  The  lateral  mass 
transport  is  also  calculated  for  the  layer  and  is  averaged 
over  time  periods.  Vertical  mass  transport  is  inferred  from 
the  vertical  velocity  (omega)  fields.  Time  periods  refer  to 
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the  12-h  tine  increment  between  synoptic  times.  For  exam¬ 
ple,  period  1206  refers  to  0000-1200  GMT.  Hadiai  distances 
of  6  and  10  degrees  of  latitude  specify  the  horizontal 
budget  area  and  correspond  to  360  and  600  n  mi  radial  dis¬ 
tances  from  storm  center.  These  particular  radii  were  cho¬ 
sen  after  inspection  of  radial  sections  depicting  the 
storm's  major  circulation  features.  The  inner  radius  (6) 
incorporates  the  zone  of  maximum  cyclone  circulation  associ¬ 
ated  with  the  strong  two-layer  inf low/cutf low  regime  (subse¬ 
quently  addressed).  The  outer  radius  (10),  in  contrast, 
shows  tha  weaker  storm  circulation  near  its  periphery. 

B.  POTENTIAL  TEMPERATURE  AND  STATIC  STABILITY  ANALYSIS 

The  steady  flow  of  polar  continental  air  from  the  Sibe¬ 
rian  mainland  over  the  warm  waters  of  the  western  Pacific 
was  described  in  Chapter  4.  Destabilization  of  the  boundary 
layer  in  this  manner  is  a  critical  factor  for  oceanic  cyclo¬ 
genesis,  as  demonstrated  by  Roman  (1981)  and  Sandgathe 
(1981).  The  amount  of  destabilization  is  contingent  upon 
tha  latent  heat  release  in  the  layer  and  the  sensible  heat 
flux  from  the  ccean.  These  quantities  are  dependent  tc  a 
large  degree  on  the  air-sea  temperature  difference,  and  the 
magnitude  cf  the  lapse  rate  through  the  lower  troposphere. 
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Roman  (1981)  found  stability  decreases  of  nearly  2  K  /  100 
mb  fully  12  h  before  the  onset  of  maximum  deepening  using  a 
budget  volume  analysis.  Numerical  model  simulations  by 
Sandgathe  (1981)  showed  a  similar  trend,  with  magnitudes 
roughly  half  that  observed  by  Roman  (198  1)  for  a  western 
Atlantic  storm.  An  analysis  of  the  destabilization  of  the 
lower  troposphere  in  this  case  will  be  described  using  the 
potential  temperature  vertical-time  sections. 

Inspection  of  the  time  sections  of  area-averaged 
potential  temperature  for  radii  6  and  10  (Fig.  27a  and  b) 
shows  similar  interesting  features.  Notice  that  potential 
temperature  values  in  this  figure  are  presented  for  12-h 
synoptic  times.  Lower  tropospheric  (1000-750  mb)  and  upper 
tropospheric  (45Q-175  mb)  potential  temperature  changes 
imply  warming  during  the  early  time  periods.  In  the  mid- 
troposphere  (750-450  mb) ,  potential  temperature  changes  lit¬ 
tle  in  the  early  time  periods.  Of  interest  is  the  apparent 
lack  of  a  well-defined  trcpopause  (indicated  by  a  dense 
vertical  packing  of  isentrcpes)  during  the  early  time  peri¬ 
ods.  Possible  factors  accounting  for  this  observation 
include  the  iirpact  of  jet  streak  propagation  (described  in 
Chapter  4)  into  the  budget  volume,  or  the  characteristic 
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upward  slops  cf  the  tropopause  toward  the  tropics,  or  a  '■ci- 
bination  of  both.  At  later  time  periods,  the  entire  budget 
volune  is  seen  to  cool,  with  the  most  significant  decrease 
occurring  below  350  mb,  as  the  cyclone  takes  on  familiar 
extr atrop ical  ccld-core  characteristics.  The  tropopause 
becomes  well  established  by  0000  GMT  15  January,  which 
coincides  with  the  movement  of  the  budget  volume  cut  of  the 
region  of  the  taroclinic  2cne  and  into  the  colder  polar  air. 
Several  substantial  potential  temperature  fluctuations  exist 
during  the  1800  GMT  times  (not  shown)  and  are  addressed 
la  ter. 

To  convert  potential  temperature  changes  to  static  sta¬ 
bility  features,  a  procedure  similar  to  that  employed  by 
Sandgathe  (1981)  will  be  used.  He  defines  the  static  sta¬ 
bility  of  the  layer  as  directly  proportional  to  the  specific 
volume  of  the  layer  and  tc  the  change  in  potential  tempera¬ 
ture  between  pressure  surfaces  (de/dP),  and  inversely  pro¬ 
portional  to  the  average  potential  temperature  of  the  layer. 
For  periods  cf  interest  here,  the  static  stability  of  the 
1000-500  mb  layer  will  be  explored.  The  specific  volume  of 
the  layer  is  assumed  constant  with  time  (variations  are  less 
than  3%)  .  The  potential  temperature  fields  (Fig.  27)  imply 
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rather  uniform  average  temperature  structure  for  the  lover 
troposphere  through  time.  Thus,  the  static  stability  varia¬ 
tions  are  predoaina ntly  a  function  cf  the  change  in  the 
lapse  rate  with  time. 

The  static  stability  trends  for  inner  and  outer  radii 
are  shown  in  Fig.  28a.  Destabilization  of  the  layer  occurs 
for  tcth  radii  until  1200  GMT  14  January.  The  increase  in 
stability  after  1200  GMT  14  January  is  a  reflection  of  cold 
advection  in  the  lower  troposphere.  The  narked  decrease  in 
stability  during  the  early  periods  (Fig.  28a)  is  associated 
with  an  increase  in  the  surface  layer  potential  temperature. 
Thus,  surface  instability  promotes  the  generation  of  organ¬ 
ized  convection  in  the  stcrm  region,  and  serves  to  foster 
the  cyclonic  circulation  tendency.  k  similar  trend  was 
observed  by  Sandgathe  (1981)  in  his  numerical  simulations 
(Fig.  28b).  Direct  comparison  of  observed  and  simulated 
stability  magnitudes  is  not  feasible  due  to  the  significant 
differences  in  the  horizontal  scales  over  whica  the  data  was 
averaged.  It  can  be  seen,  however,  that  the  trends  in  the 
observat icnal  data  match  the  simulated  data  very  closely, 
and  show  the  characteristics  of  one  of  the  main  ingredients 
that  is  favorable  for  oceanic  cy clo genesis . 
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Cns  interesting  feature  of  Fig.  27a  deserves  special 
scrutiny.  The  dramatic  fluctuations  in  potential  tempera¬ 
ture  (8  K)  between  1200  GBT  14  January  and  0000  GMT  15  Janu¬ 
ary  between  8C0  and  500  mb  is  remarkable  for  such  a  small 
time  interval.  In  general,  the  atmosphere  requires  much 
longer  periods  fcr  such  large  fluctuations.  Inspection  of 
85  0-500  mb  potential  temperature  fields  for  both  radii  for 
times  immediately  preceding  and  subsequent  to  1200  GMT  shows 
significant  warming  evident  in  the  1800  GMT  field.  In  addi¬ 
tion,  potential  temperature  biases  of  smaller  magnitudes 
were  evident  at  the  1800  and  0600  GMT  times  on  other  days 
(Fig.  27a  and  b) .  The  nature  of  these  (1800  and  0600  GMT) 
potential  temperature  fluctuations  suggests  a  problem  with 
the  FGGE  data  assimilation  schemes,  and  its  impact  on  this 
study  is  specifically  addressed  in  Appendix  A. 

C.  HCRIZCSTAI  MASS  FLUX 

The  system  horizontal  mass  transport  analysis  consists 
of  vertical  time  sections  for  radii  6  and  ’0  (Fig.  29a  and 
b)  in  terms  of  g/sec  for  the  100  mb  layer.  Flux  values  are 
corrected  for  mean  cyclone  movement  and  residual  imbalances 
(see  Chapter  3).  Time  periods  refer  to  the  12-h  period 
between  synoptic  times.  For  example,  period  1206  refers  to 
0000-1200  GMT  12  January. 
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Inspection  of  Pig.  29a  for  radius  six  shows  a  two-layer 
regime  consisting  of  strong  low-level  mass  influx  associated 
with  surface  layer  convergence  below  a  layer  of  mass  outflow 
associated  with  middle  and  upper-level  divergence.  The 
strength  cf  the  mass  flew  is  seen  to  increase  in  both 
regimes  as  the  storm  develops,  with  a  distinct  inflow  maxi¬ 
mum  near  925  mb  and  an  outflow  maximum  near  350  mb.  The 
inflow  maximum  (3.27  X  10  *  *  13)  at  13  18  GUT  is  roughly 
twice  the  cutflcw  maximum  (1.56  X  10  *  *  13).  Although  the 
low-level  inflow  is  intense,  it  is  restricted  to  a  rela¬ 
tively  narrow  layer.  The  deeper,  though  weaker,  upper-level 
outflow  causes  a  net  vertically  integrated  mass  loss  within 
the  column.  The  net  horizontal  mass  outflow  in  the  vertical 
columns  (Fig.  29a)  is  reflected  in  the  surface  pressure 
falls  (Fig.  2b)  during  the  developing  periods. 

The  temporal  variation  in  the  level  of  non-divergence 
(LND)  and  the  associated  inflow  maxima  (Fig.  29a)  between 
1306  and  1406  3  MT  indicates  vertical  growth  of  the 
convergent  layer  as  intense  cyclogenesis  occurs.  3y  1306 
GMT,  when  closed  isobars  become  evident  at  the  surface,  the 
convergent  layer  begins  tc  grow  as  the  surface  circulation 
is  established.  The  growth  of  the  convergent  layer  during 
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cyclcgenesis  is  consistent  with  results  found  by  Soman 
(198  1).  The  radius  6  LND  reaches  maximum  height  (at  1406 
GMT)  soon  after  the  occurrence  of  the  maximum  surface  layer 
transport  and  just  before  the  period  of  lowest  surface  pres¬ 
sure  (1500  GMT),  and  then  exhibits  a  gradual  reduction  with 
time  as  the  stern  begins  to  dissipate. 

Comparable  features  exist  within  radius  10  (Fig.  29b) 
but  interesting  differences  exist.  The  mass  buildup  in  the 
outer  volume  (1318  GMT)  shows  a  greater  relative  increase, 
which  indicates  differential  transport  between  the  two 
radii,  since  the  inner  volume  increases  are  incorporated  in 
the  cuter  volume  tendency.  The  differential  mass  transport 
car.  be  explained  by  the  earlier  influence  of  the  migratory 
jet  streak  (discussed  in  Chapter  4)  as  it  propagates  first 
into  the  cuter  budget  volume.  The  decrease  in  mass 
transport  after  1318  GMT  reflects  the  decrease  in  jet  streak 
influence  as  it  propagates  through  the  inner  volume  and  dis¬ 
sipates  (see  Chapter  4).  These  patterns  substantiate  the 
intensity  cf  associated  "weather"  observed  by  stations  hun¬ 
dreds  of  n  mi  from  the  storm  center  and  clearly  illustrate 
the  immense  scale  of  these  meteorological  events.  It  should 
also  be  noticed  that  the  increase  in  low-level  circulation 
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during  the  last  few  time  periods  at  radius  10  (Fig.  29b)  is 
attributed  tc  the  incorporation  of  the  developing 
circulation  to  the  south  (see  Chapter  4). 

Several  interesting  features  of  the  development  stage  of 
this  cyclone  can  be  observed  by  inspection  of  the  radial- 
vertical  sections  of  the  corrected  lateral  mass  flux.  A 
sequence  cf  fcur  times  was  chosen  (0000  and  1200  GMT  14  Jan¬ 
uary  and  0000  and  1200  GMT  15  January)  illustrated  in  Figs. 
30  and  31  respectively. 

At  the  outer  radii  in  Fig.  30a,  the  mass  transport  is  a 
maximum  near  radius  11.  The  LND  is  relatively  constant  in 
height  through  all  radii.  A  broad  region  of  strong  upper 
level  divergence  is  also  evident  near  the  outer  radius. 
Notice  the  significant  change  twelve  hours  later  (Fig. 
30b).  An  additional  layer  of  surface  convergence  near  775 
mb  approaches  the  storm  center  from  the  outer  radii.  This 
feature  suggests  that  a  major  significant  change  in  the 
circulation  has  occurred  during  the  past  twelve  hours.  Nor¬ 
mal  cyclone  development  might  explain  a  portion  cf  the 
increase,  but  the  explosive  development  indicated  during 
this  time  frcm  the  rapid  pressure  falls  suggests  seme 
additional  forcing  present  in  the  budget  volume. 
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In  the  subsequent  1 2  —  h  period  (Fig.  31a),  the  surface 
transport  maximum  decreases  and  oo7es  well  inside  radius  6. 
Divergence  aloft  is  well  established  and  the  positive 
transport  inside  radius  3  between  400-125  mb  suggests  net 
subsidence  and  adiabatic  warming  near  the  cyclone  core, 
which  arrests  at  least  to  some  small  degree  the  cold-ccre 
formation  process.  Transport  values  markedly  decrease  12  h 
later  (Fig.  31b)  and  suggest  cyclone  dissipation. 

These  circulation  features  may  be  compared  with  major 
characteristics  cf  tropical  storm  intensification  (e.g.  Hol¬ 
land  and  Merrill,  1  983).  Explosive  development  in  both 
cases  appears  to  be  char acteriz ed  by  large  circulation  and 
mass  transport  increases  at  the  outer  radii,  with  subsequent 
inner  radial  development  after  organized  transport  has  moved 
near  the  inner  region  with  deep  convection.  The  inner 
radial  development  is  also  largely  dependent  on  the 
formation  of  a  broad  and  strong  outflow  regime  in  the  upper 
troposphere.  As  these  significant  features  form,  organized 
convection  dominates  the  inner  radii  and  triggers  the  rapid 
pressure  falls  due  to  large  latent  heat  inputs  into  the 
lower  troposphere.  Large  pressure  falls  set  the  stage  for 
increased  winds,  circulation  etc..  The  storm  rapidly 
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develops  with  warming  near  the  central  core  until  the 
surface  system  becomes  detached  from  the  upper-level  venti¬ 
lation  and  divergence  source.  The  significant  observed  dif¬ 
ference  in  the  formation  processes  lies  in  the  expanded 
radial  scale  of  the  extratropi cal  system  compared  to  the 
mesoscale  structure  of  the  tropical  system.  Otherwise, 
upper-level  ventilation  features,  inner  radius  organized 
deep  convection,  weak  stability  of  the  surface  layer  and 
large  outer  region  transports  in  the  early  development  peri¬ 
ods  suggest  reasonable  similarities.  The  conclusions  cited 
above  are  only  suggestions  from  the  available  FGGE  fields 
and  will  not  be  explored  in  detail,  as  these  aspects  are  net 
considered  a  primary  objective  of  this  study. 

D.  VERTICAL  MASS  TRANSPORT 

The  storm  vertical  mass  transport  analysis  consists  cf 
vertical  time  sections  of  cmega  fields  (Figs.  32,  33  and  34) 
for  radius  6  and  10.  Three  distinct  omega  fields  are  calcu¬ 
lated  by  slightly  different  methods  and  the  differences  are 
examined.  The  omega  fields  discussed  below  represent  the 
vertical  mass  transport  and  necessarily  reflect  periods  of 
intense  hcrizcntal  mass  transport  through  continuity  (see 
Chapter  3) . 
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The  differences  between  the  three  omega  fields  are  shown 
clearly  in  Figs.  32,  33  and  3U.  The  budget  calculated 
(QLD)  omega  fields  (Fig.  32)  are  area-averaged  values 
derived  from  the  lateral  transport  through  mass  continuity. 
The  FGGE  Level  Ill-b  vertical  motion  estimates  (Fig.  33)  are 
derived  as  point  estimates  without  diabatic  internal  adjust¬ 
ments.  Omega  fields  calculated  at  grid  points  using  the 
kinematic  method,  which  implicitly  includes  small-scale 
diabatic  effects,  are  shown  in  (Fig.  34). 

The  difference  in  magnitude  and  level  of  the  vertical 
motion  maxima  between  the  three  vertical  motion  fields  is 
interesting.  The  QLD  (Fig.  32)  and  FGGE  (Fig.  33)  vertical 
velocity  maxima  match  fairly  closely  in  time  of  occurrence 
and  magnitude,  but  differ  considerably  in  level.  Compensa¬ 
tion  for  the  residual  imbalances  in  the  budget  calculations 
(discussed  in  Chapter  3)  are  applied  layer  by  layer  (highest 
layers  are  weighted  more)  in  the  QLD  omega  field  derivation. 
Thus,  the  QLD  omega  magnitudes  are  comparable  but  the  level 
of  the  maxima  is  biased  toward  the  higher  layers  (Fig.  32). 

The  kinematic  vertical  velocities  differ  both  in  magni¬ 
tude  and  level.  This  is  not  surprising  considering  the  val¬ 
ues  include  small-scale  diabatic  effects.  Paegle  (1983) 
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demonstrated  similar  results  and  found  that  the  FGGE  derived 


fields  were  roughly  1/3  the  magnitude  of  the  kinematically 
derived  fields,  which  is  consistent  with  a  comparison  of 
Figs.  33  and  34  (the  ratio  is  closer  to  2/3  here).  For  the 
purpose  of  this  study,  the  kinematic  omega  fields  are 
considered  to  test  represent  the  actual  vertical  motion. 

h  comparison  of  the  vertical  velocity  fields  with  the 
horizontal  mass  transport  shows  consistent  features.  The 
vertical  velocity  maxima  (Fig.  34a)  correlate  well  tempo¬ 
rally  with  the  surface  convergence  and  upper  level 
divergence  maxima  (Fig.  29a),  and  with  periods  of  rapid 
surface  deepening  (Fig.  2t) .  In  addition,  the  more  intense 
vertical  motion  at  radius  6  (Figs.  32,  33  and  34) ,  which  is 
roughly  1.5  times  greater  than  at  radius  10,  agrees  well 
with  the  satellite  discussion  which  depicted  extreme  amounts 
of  deep  convection  in  the  inner  radius  (see  Chapter  4) . 

E.  CCNCLU  SICKS 

In  summary,  stability  and  QLD  mass  budget  analyses  have 
revealed  the  following  insights: 

•  Surface  stability  decreases  during  the  early  time  peri¬ 
ods  favor  organized  convection; 

•  The  stability  analysis  shows  destabilization  trends 
consistent  with  numerical  model  simulations  (Sandgathe, 
1981); 
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•  Evolution  of  cold  core  character istics  is  evident 
throughout  tost  of  the  troposphere  after  the  surface 
circulation  intensifies; 


•  Maxima  in  surface  layer  convergence  occur  during  peri¬ 
ods  of  intense  vertical  motion  and  strong  upper-level 
divergence,  coincident  with  observed  maximum  cyclone 
deepening  rates  (Fig.  2b); 


•  The  bread  divergent  layer  aloft  contributes  to  central 
pressure  falls  even  though  the  maxima  are  considerably 
smaller  than  the  associated  surface  layer  convergent 
maxima,  indicating  overall  net  mass  outflow  during 
periods  of  storm  intensification; 


•  The  LND  varies  according  to  the  strength  of  the  surface 
convergent  layer,  and  rises  during  the  period  of  explo¬ 
sive  development; 


•  Vertical  velocity  fields  employing  Kinematic  calcula¬ 
tion  methods  represent  well  the  LMD  and  vertical  mass 
transport  which  is  a  vital  link  in  the  lifeblood  of 
storm  dynamics. 


•  Vertical-radial  cross  sections  of  mass  transport  during 
extratrcpjca}  cyclogenes4s  suggest  similar  development 
characteristics  to  tropical  storm  rapid  intensifica¬ 
tion. 


The  next  chapter  seeks  to  explore  the  circulation 


tendencies  as  delineated  by  the  absolute  vorticity  buildup 


within  the  budget  volume  as  the  storm  evolves.  Horizontal 


and  vertical  transport,  and  sources  and  sinks 
vorticity  will  be  explored  guantitati vel y 
delineate  storm  dynamical  features. 


of  absolute 
to  further 
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VI.  VOETICITY  BOD GET  ANALYSIS 

A.  GENERAL 

The  absolute  vorticity  budget  analyses  are  presented  by 
vertical-time  sections  and  by  selected  vertical  profiles. 
Time  pericds  refer  to  the  12- h  span  between  synoptic  times. 
For  example,  period  1206  refers  to  0000-1200  GMT  12  January. 
In  figures  not  referencing  time  periods,  times  refer  to  the 
specific  time  (GMT)  .  The  775  mb  layer  refers  to  the  layer 
between  the  650  and  700  mb.  The  sign  convention  is  positive 
for  processes  producing  vcrticity  increases  in  the  budget 
volume  (sources),  and  negative  for  vorticity  decreases 
(sinks)  . 

The  main  feature  of  this  section  includes  a  discussion 
of  the  absolute  vorticity  time  tendency  as  a  measure  of 
stcrm  development.  Specifically,  an  investigation  is  con¬ 
ducted  intc  the  role  of  horizontal  and  vertical  absolute 
vorticity  transport,  advection,  generation  and  dissipation 
in  storm  dynamics.  Horizontal  transport  terms  are  separated 
into  mean  and  eddy  mode  components  to  gain  better  insight 
into  their  specific  contribution s  toward  storm  development. 
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Vertical  transport  terns  reflect  vertical  divergence  and 
advection  frca  caega  estinates  using  the  kinematic  method. 
Referral  to  section  B  of  Table  V  is  recommended  as  specific 
budget  terms  are  identified  and  discussed. 

B.  TIME  SECTIONS  ,  TIME  TENDENCIES  AND  SELECTED  PROFILES 

The  dramatic  buildup  of  absolute  vorticity  in  the  budget 
volume  at  both  radii  over  the  2-day  span  of  intense 
development  is  illustrated  in  Fig.  35.  At  0 000  GMT  13  Janu¬ 
ary  (1300  GMT) ,  the  relative  cyclonic  vorticity  in  the  inner 
volume  through  the  600  mb  level  is  not  much  different  from 
the  ccriolis  farameter  (7.82  X  10**-5).  Notice  the  presence 
of  anticyclonic  relative  vorticity  aloft,  which  is  associ¬ 
ated  with  the  weak  and  dissipating  upper  level  ridge.  The 
dramatic  increase  in  absolute  vorticity  by  0000  GMT  15  Janu¬ 
ary  in  the  inner  radius  occurs  from  two  sources.  The  rapid 
movement  to  the  northeast  over  the  period  has  resulted  in  a 
larger  value  of  planetary  vorticity  (f  =  11.56  Xl0**-5), 
while  the  increased  storm  circulation  accounts  for  the 
additional  buildup.  Although  the  surface  circulation  is 
strong,  the  iraxiaua  absolute  vorticity  exists  in  the  middle 
troposphere.  Similar  changes  occur  at  radius  10,  although 
the  magnitude  of  the  absolute  vorticity  increase  is  smaller. 
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and  there  is  a  uniform  increase  through  the  troposphere 
below  the  jet  core  level.  Specific  factors  contributing  to 
the  vorticity  changes  at  the  different  vertical  levels  will 
be  discussed  as  the  specific  budget  terms  are  analyzed. 

A  different  perspective  of  the  absolute  vorticity 
buildup  in  the  budget  volume  is  presented  in  Fig.  36.  Here, 
vertical  time  sections  of  vorticity  for  both  radii  illus¬ 
trate  several  key  features.  At  radius  6  (Fig.  36a),  a  vor¬ 
ticity  minimus  can  be  seen  at  200  mb  during  the  early 
period.  The  increase  with  time  at  lower  levels  is  a  reflec¬ 
tion  of  the  nascent  cyclone  vertical  structure  which  is  com¬ 
posed  of  anticyclonic  vcrticity  aloft  over  the  low-level 
developing  system.  After  1212  GMT,  absolute  vorticity 
increases  throughout  the  whole  troposphere  until  1500  GMT. 
Two  significant  features  are  evident  during  these  rimes.  A 
large  buildup  occurs  in  the  upper  troposphere  (vicinity  225 
mb)  after  1300  GMT,  which  causes  a  remporary  reduction  in 
the  vertical  vorticity  gradient.  A  distinct  lower  tropos¬ 
pheric  maximum  is  evident  at  1500  GMT,  which  occurs  just 
after  the  period  of  intense  circulation  increase,  and  is 
substantiated  by  the  mass  transport  maxima  shown  in  the  mass 
budget  analysis  (Fig.  29a)  and  the  associated  central 
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surface  pressure  falls  illustrated  in  Fig.  2a  of  the  synop¬ 
tic  discussion.  This  large  maximum  in  vorticity  is  coinci¬ 
dent  with  the  period  of  lowest  central  pressure  and  maximum 
surface  layer  circulation.  Vorticity  values  are  seen  to 
decrease  in  the  middle  and  lower  troposphere  for  24  h  after 
1500  GMT,  but  at  a  significantly  slower  rata  than  the 
buildup.  Values  increase  again  through  the  whole 
troposphere  after  1600  GMT  due  to  the  inclusion  of  the 
developing  system  to  the  scuth. 

Several  important  differences  exist  at  radius  10  (Fig. 
36b).  vorticity  values  increase  at  a  slower  rate  and  at 
later  times  than  in  the  inner  volume.  This  result  indicates 
that  although  vorticity  values  increase  in  both  volumes  dur¬ 
ing  the  development  stage,  the  outer  volume  experiences  a 
less  significant  increase  due  to  the  greater  spatial  extent 
and  relatively  smaller  vorticity. 

Absolute  vorticity  time  tendency,  which  is  the  first 
term  in  the  vorticity  budget  equation  (see  Table  V) ,  is  pre¬ 
sented  in  Fig.  37.  This  display  is  a  time  derivative  of  the 
fields  depicted  in  Fig.  36.  Recall  that  time  periods,  when 
used  in  subsequent  figures,  refer  to  tHa  12-h  time  interval 
between  synoptic  times. 
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A  large  increase  in  tie  vorticity  tendency  in  the  inner 
volume  occurs  in  the  upper  levels  from  1218-1318  GMT.  Vor- 
ticity  increases  at  a  constant  rate  in  the  surface  layer 
during  the  development  stage.  The  general  decrease  in  ten¬ 
dency  values  through  the  whole  troposphere  after  1506  GMT  is 
consistent  with  the  onset  of  the  occlusion  phase  and  storm 
dissipa1  ion.  It  is  interesting  that  a  vorticity  tendency 
minimum  occurs  near  450  mb  during  a  period  of  large  surface 
pressure  falls  (Fig.  2b)  and  maximum  mass  transport  (Fig. 
29).  No  correlation  to  synoptic  causes  for  this  tendency 
minimum  is  apparent  from  detailed  inspection  of  the  500  mb 
height  and  vorticity  fields  (Figs.  4b,  7b,  10b).  One  inter¬ 
pretation  is  tnat  a  subsynoptic  scale  system  in  the  wind 
fields  affects  the  vorticity  values  at  these  levels  and 
times.  An  alternative  cause  could  be  United  to  insufficient 
data  input  into  the  1312  or  1400  GMT  analyses. 

Similar  significant  features  exist  at  radius  10  (Fig. 
37b).  The  nature  of  the  buildup  of  absolute  vorticity 
suggests  inherently  similar  mechanisms  are  responsible  for 
the  absolute  vcrticity  buildup  between  inner  and  outer  vol¬ 
umes.  Mechanisms  such  as  the  jet  core  influence,  mid-trc- 
pospheric  self-acplif icaticn,  and  horizontal  and  vertical 
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transport  are  possibilities  explored  in  the  following 
sections. 

C.  LATERAL  TRANSPORT 

The  lateral  transport  of  absolute  vorticity  is  one  of 
the  three  most  significant  forcing  terms  in  the  budget  equa¬ 
tion  (Table  V)  .  Transport  into  the  budget  volume  arises 
from  the  mean  convergent  flow  (mean  mode)  and  from 
asymmetries  in  the  cyclonic  flow  (eddy  mode) . 

A  typical  view  of  the  lateral  vorticity  transport  occur¬ 
ring  during  storm  development  is  presented  in  Fig.  38.  The 
curves  are  valid  during  the  maximum  development  period  at 
radius  six.  The  mean  mode  and  total  lateral  transport 
curves  are  similar  and  shew  similar  trends.  The  mean  mode 
is  seen  to  dominate  the  total  transport  at  levels  below  800 
mb.  In  the  upper  levels,  the  total  lateral  transport  is 
small,  but  is  the  residual  between  two  strong  oppositely 
contributing  processes.  These  features  constitute  signifi¬ 
cant  differences  from  those  observed  by  Conant  (1982)  in  the 
President's  Day  storm,  where  he  observed  that  the  total  and 
eddy  mode  lateral  transport  curves  were  the  more  strongly 
correlated.  Although  the  curves  presented  in  that  study 
represented  an  earlier  phase  in  the  storm  development,  the 
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general  features  indicate  that  different  mechanisms  serve  to 
increase  inner  volume  vorticity  levels  within  each  storm 
budget  volume.  A  look  at  the  mean  and  eddy  mode  lateral 
transport  vertical-time  sections  gives  insights  into 
processes  responsible  for  these  differences. 

Another  interesting  feature  is  the  presence  of  a  measur¬ 
able  eddy  mode  contribution  in  the  lower  troposphere,  which 
is  an  indication  that  processes  other  than  mean  convergence 
are  acting  to  increase  the  vorticity  and  storm  circulation 
within  the  six  degree  volume.  These  results  should  be 
expected,  since  processes  such  as  intense  convection  in  pre¬ 
ferred  quadrants  of  the  storm,  frontal  asymmetries  and 
elliptical  vertex  structure  may  cause  asymmetric  flew  about 
and  into  the  cyclone.  Thus,  some  eddy  mode  contribution  is 
evident  through  the  whole  troposphere.  Mean  and  eddy  mode 
contributions  to  the  total  vorticity  lateral  transport  are 
discussed  belcw. 

The  xotal  lateral  vorticity  transport  rime  sections  for 
the  inner  and  outer  budget  volumes  are  presented  in  Fig. 
39.  In  the  inner  volume  (Fig.  39a),  inward  transport  exists 
below  775  mb  through  all  time  periods,  with  alternating 
inward  and  outward  transport  evident  in  the  upper 
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troposphere  (near  350  mb)  .  Notice  the  mid-tropospheric  (600 
mb)  outward  transport  (after  1506  GMT)  as  broad  divergence 
is  established.  The  outward  transport  maximum  at  225  mb 
suggests  strong  divergence  aloft  as  the  jet  stream  dominates 
the  upper  troposphere  flew.  The  surface  layer  regime  is 
particularly  noteworthy  in  that  especially  deep  inward  vor- 
ticity  transport  occurs  during  the  deepening  phase 
(1318-1506  GMT)  through  450  mb  and  contains  a  maximum  (19.2 
X  10**-10/sec**2)  at  925  mb.  This  signature  is  entirely 
different  from  that  observed  by  Conant  (1982)  in  the  Presi¬ 
dent's  Day  storm.  He  shewed  a  maximum  at  220  mb  with  rela¬ 
tively  weak  surface  transport.  The  emphasis  on  the  level  of 
maximum  hcrizcttal  lateral  transport  serves  to  focus  on  the 
differences  between  the  underlying  processes  responsible  for 
vorticity  increases  within  the  budgets  of  both  systems.  In 
this  case,  the  lateral  transport  surface  maximum  occurs  in 
temporal  coincidence  with  the  largest  central  pressure  falls 
and  maximum  mass  transport,  as  discussed  in  previous 
chapters. 

The  tctal  lateral  transport  at  radius  10  (Fig.  39b) 
exhibits  similar  inward  transport  in  the  surface  layer 
although  the  magnitude  is  considerably  less  and  the  maximum 
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occurs  a  time  period  earlier.  Inward  vorticity  transport  is 
evident  through  the  rest  cf  the  troposphere  (except  above 
450  mb  after  1606  GMT) ,  with  a  secondary  maximum  observed 
between  450  and  175  mb  during  the  1218-1318  GMT  periods. 
The  time  and  level  of  occurrence  of  this  maximum  strongly 
suggests  vcrticity  import  from  jet  streak  interaction.  Sub¬ 
sequent  low-level  inward  transport  occurs  from  circulation 
increases  in  response  tc  this  upper-level  forcing.  The 
earlier  occurrence  of  both  upper  and  lower-level  maxima  in 
the  cuter  volume  suggests  earlier  forcing  from  jet  streak 
interaction  as  it  enters  the  outer  volume  first  and  subse¬ 
quently  migrates  toward  the  inner  volume.  The  increase  in 
the  surface  layer  transport  at  1618  GMT  is  again  due  to  the 
developing  system  to  the  scuth. 

C.  MEAN  MCDE  LAT2BAL  TRANSPORT 

As  described  in  Chapter  3  (Table  V)  ,  the  absolute  vcr¬ 
ticity  transpcrt  can  be  partitioned  into  mean  and  eddy 
modes.  The  purpose  of  this  partitioning  serves  to  isolate 
contributions  by  the  mean  flow  (mean  cyclone  convergence  and 
divergence),  and  the  perturbation  flow  (the  covariance  of 
vorticity  and  normal  wind  component  deviations) . 
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Specifically,  the  mean  cyclonic  low-level  inflow  and  upper- 
level  outflow  cf  vorticity  is  represented  by  the  mean  mode 
and  reflects  necessarily  the  convergent/divergent  nature  of 
the  lcver/upper  layers. 

At  radius  six  (Fig.  40a) ,  the  strong  dual-centered 
transport  maximum,  the  upper-level  outflow  maxima  and  the 
level  of  transition  between  inward  and  outward  transport 
agree  well  both  temporally  and  spatially  with  the  mass 
budget  transport  features  discussed  in  Chapter  5.  Periods 
of  maximum  surface  cyclone  deepening  coincide  directly  with 
the  vorticity  transport  maxima  discussed  above.  Differences 
in  mean  mode  transport  between  the  President's  Day  storm  and 
the  present  case  arise  due  to  the  significantly  weaker 
surface-layer  lateral  vorticity  transport  in  the  President's 
Day  storm,  whereas  the  upper  layer  transports  are  nearly  the 
same.  Thus,  the  main  difference  in  the  growth  and  vitality 
in  this  storm  is  the  significantly  higher  mean  flow 
transport  of  vorticity  in  the  lower  layers. 

At  radius  10  (Fig.  40b),  absolute  vorticity  lateral 
transport  values  are  relatively  weaker  than  at  radius  6.  A 
single  inward  transport  maximum  in  the  surface  layer  and  a 
coincident  upper-level  outward  maximum  are  evident.  These 
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maxima  occur  earlier,  again,  in  response  to  earlier  dynamic 
forcing  in  the  outer  vclume  by  the  upper-level  eddy 
structure  changes  (subsequently  addressed) .  The  effects  of 
the  new  lew  tc  the  south  are  again  seen  after  1518  GMT. 

E.  EDDY  MODE  LATERAL  TRANSPORT 

The  eddy  mode  absolute  vorticity  lateral  transport  gives 
insight  into  the  upper  level  vorticity  increases  which  ulti¬ 
mately  affect  the  lew  level  circulation  tendencies.  As 
discussed  above,  the  eddy  mode  component  represents  the  hor¬ 
izontal  transport  due  to  asymmetries  in  the  cyclonic  flow 
induced  by  wind  field  correlated  vorticity  deviations. 

Inspection  of  the  eddy  mode  absolute  vorticity  transport 
for  radius  6  (Fig.  4  1a)  shows  considerable  vorticity  import 
in  the  layer  between  750  and  125  mb  after  1206  GMT.  within 
this  broad  area  of  positive  lateral  vorticity  transport  are 
two  significant  time  periods.  Between  1206-1306  GMT,  a  max¬ 
imum  of  vorticity  transport  is  seen  to  exist  in  the  450  275 
mb  layer.  Inspection  of  the  300  and  400  mb  vorticity  and 
winds  for  1200  and  1212  GMT  shows  weak  PVA  into  the  budget 
volume,  and  this  is  substantiated  by  the  coincident  maximum 
in  the  calculated  advection  (Pig.  43a)  discussed  subse¬ 
quently.  These  facts  suggest  vorticity  import  in  the  form 
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of  the  short  nave  that  becomes  established  over  southern 


Japan  during  this  time  frame  (discussed  earlier  in  the 
synoptic  discussion.  Fig.  4b).  The  effect  of  this  perturba¬ 
tion  diminishes  with  time  as  the  cyclone  accelerates  to  the 
northeast  away  from  the  slower  moving  upper-level  trough 
after  1300  GMT,  and  thus,  the  eddy  mode  transport  also 
diminishes.  Vorticity  increases  are  again  evident  after 
1306  GMT  when  the  jet  streak  has  become  well  established  in 
the  budget  volume  (Fig.  7a).  The  influence  of  the  jet 
streak  continues  until  1518  GMT  when  the  jet  core  becomes 
detachad  from  the  budget  volume  as  the  surface  center 
migrates  steadily  northward. 

Significant  eddy  mode  contributions  are  restricted  to 
the  upper  levels,  and  are  only  roughly  1/2  those  observed  by 
Conant  (1982)  in  the  President’s  Day  storm.  This  is  an 
extremely  valuable  result,  as  Ccnant  showed  that  the  migra¬ 
tory  short  wave  played  a  major  role  in  that  storm 
development.  In  this  study,  the  migratory  short  wave  is 
shown  to  be  of  minimal  influence,  and  thus  vorticity 
increases  in  the  upper  levels  occur  primarily  from  cyclonic 
shear  vorticity  input  into  the  budget  volume  as  the  jet  core 
becomes  established.  While  the  eddy  mode  contributes 
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considerably  to  the  initial  storm  development,  this  affect 
is  net  nearly  as  dominant  during  the  explosive  stages  as  in 
the  President's  Eay  storm.  In  this  case  then,  it  is  logical 
to  conclude  that  explosive  development  occurs  due  to  upper- 
level  jet  forcing  acting  in  concert  with  the  surface 
development  alcng  the  barcclinic  zone.  The  secondary  eddy 
mode  maximum  near  925  mb  at  1418  GMT  is  presumed  to  be  due 
to  asymmetries  ir.  the  mean  flow  in  specific  quadrants  of  the 
storm  circulation. 

Similar  features  exist  at  radius  10  (Fig.  41b).  As  in 
the  inner  volume,  significant  contributions  occur  at  the  jet 
core  level  with  maxima  cf  similar  magnitude  occurring  at 
coincident  times.  A  weak  maximum  exists  from  1406-1506  GMT 
in  the  surface  layer,  and  as  in  the  6-degree  volume, 
represents  small  asymmetries  in  the  mean  cyclonic  flow. 

Alternatively,  the  lateral  transport  may  be  expressed  in 
terms  of  lateral  divergence  and  advection  components  through 
the  use  of  vector  identitites  (see  Table  V).  These  terms 
provide  a  more  traditional  view  directly  in  terms  of  the 
vorticity  equation.  The  divergence  component  of  the  parti¬ 
tioned  lateral  vorticity  transport  (Fig.  42a),  is  composed 
of  a  two-layer  regime  with  positive  vorticity  production  due 
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to  convergence  in  the  surface  layers  and  negative  vorticity 
production  with  divergence  aloft.  .  These  features  match 
closely  the  associated  mass  transport  features  discussed  in 
Chapter  5.  The  temporal  and  spatial  occurrence  of  the 
inflow  and  outflow  maxima  agree  well  with  similar  features 
of  the  mean  mode  discussion,  and  with  the  general 
circulation  increases  and  maximum  pressure  falls. 

The  most  interesting  comparison  occurs  between  the  eddy 
mode  and  hcrizcntal  advection  transports  (Figs.  41a  and  b 
and  43a  and  fc)  .  The  upper  tropospheric  values  match  well 
both  temporally  and  spatially,  as  well  as  in  magnitude. 
However,  the  advection  term  in  the  lower  troposphere  shows  a 
substantial  negative  vorticity  advection  (NVA)  into  the  vol¬ 
ume  below  775  mb,  which  is  quite  different  from  the  negative 
and  positive  eddy  vorticity  transports.  In-depth  examina¬ 
tion  of  the  vcrticity  and  wind  fields  in  the  inner  volume 
during  a  period  of  maximum  development  shows  an  asymmetric 
east-west  oriented  vorticity  field  coupled  to  a  strong 
cyclcnic  inflow  pattern,  which  causes  a  net  advection  of 
smaller  vorticity  values  into  the  volume  in  the  northeast 
and  southwest  guadrants.  This  net  HVA  is  further  enhanced 
by  the  40-kt  translation  speed  of  the  surface  center  to  the 
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northeast  during  this  period.  The  change  to  a  westerly, 
more  zonal  wind  component  with  height  accounts  for  the  s ion 
reversal  in  the  vorticity  values  above  775  mb.  Notice  that 
the  addition  of  the  advection  and  divergence  terms  is  very 
close  to  the  sum  of  the  mean  and  eddy  mode  lateral 
transports.  These  results  serve  to  verify  the  integrity  of 
the  transport  calculations. 

In  summary,  the  total  horizontal  absolute  vorticity 
transport  has  teen  shown  tc  be  composed  of  a  mean  mode  rep¬ 
resenting  ccntributions  from  the  mean  mass  circulation,  and 
an  eddy  mcde  representing  contributions  from  asymmetries  in 
the  flow  from  advection  and  jet  core  sources.  The  mean  mode 
transport  was  shown  to  be  very  similar  to  the  total 
transport  through  all  levels,  while  the  eddy  mode  transport 
was  shown  to  be  significant  only  above  500  mb.  The  rela¬ 
tively  smaller  magnitude  of  the  eddy  mode  transport  maxima 
compared  to  than  observed  in  the  President's  Day  storm 
suggests  this  term  play s  a  less  significant  role  in  the 
storm  evolution. 

P.  VERTICAL  REDISTRIBUTION 

This  section  investigates  the  role  of  vertical  vorticity 
transport  from  the  upper  and  lower  layers  into  the  middle 


troposphere  (as  discussed  earlier  in  connection  with  Fig. 
34).  In  this  case  the  diffluent  nature  of  the  upward 
avolving  trough  will  be  explained  in  terms  of  this  mid-level 
vorticity  increase.  The  concentration  of  vorticity  in  the 
ai d-troposphere  serves  to  increase  mid-level  circulation, 
increase  therial  advection,  and  amplify  the  trough. 
Vertical  transport  of  vorticity  was  calculated  using  the 
kinematic  omega  fields  (Chapter  5)  , 

The  affect  of  the  vertical  transport  in  the  vorticity 
budget  is  given  by  the  divergence  of  the  vertical  transport 
(Fig.  44a  and  b.)  .  This  term  results  from  the  vertical 
derivative  of  the  transport  between  levels  and  represents 
the  redistribution  of  vorticity  upward  through  the 
troposphere.  Notice  the  two-layer  regime  with  strong 
transport  from  the  surface  layer  into  the  muddle  and  upper 
troposphere.  The  vertical  transport  of  vorticity  is  a  con¬ 
sequence  of  upper-level  divergence  which  induces  lower-level 
circulation  and  vertical  xctioa  increases  throughout  the  low 
and  middle  troposphere.  The  maxima  in  both  upper  and  lower 
regimes  are  temporally  coincident  with  periods  of  maximum 
vertical  motion  (see  Chapter  5).  Values  in  the  outer  radius 
are  considerably  smaller  which  is  consistent  with  the  less 
intense  vertical  motions. 
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The  divergent  component  (Fig.  45a)  is  the  main 
contributor  tc  vertical  vorticity  redistribution.  This  term 
is  simply  the  opposite  of  the  lateral  divergence  term  (Fig. 
42a)  ,  but  differs  slightly  due  to  the  kinematic  vertical 
motion  estimates.  There  is  a  decrease  of  vorticity  below 
the  maximum  in  the  vertical  motion,  and  a  gain  above  this 
maximum.  The  reversal  in  sign  from  the  lateral  divergence 
term  (Fig.  42)  is  consistent  from  continuity,  and  close 
agreement  between  the  temporal  and  spatial  magnitudes  shows 
consistency.  Similar  features  exist  at  radius  10,  although 
the  magnitudes  are  considerably  less  as  expected. 

The  vertical  advection  portion  (Fig.  46a)  plays  a  rela¬ 
tively  minor  role  in  the  budget.  Positive  vertical  advec¬ 
tion  is  seen  tc  occur  between  775  mb  and  350  mb  from 
1318-1506  GMT.  Accumulated  values  of  surface  layer  vortic¬ 
ity  are  advected  by  the  vertical  motion  during  these  times. 
Positive  values  below  and  negative  values  above  450  mb  con¬ 
tribute  to  the  vorticity  buildup  in  the  middle  troposphere 
illustrated  in  Fig.  35.  The  relatively  minor  contribution 
of  vertical  advection  to  the  budget  results  is  consistent 
with  results  found  by  Conant  (1982). 
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In  summary,  the  vertical  transport  redistributes  the 
large  vorticity  buildup  in  the  lover  troposphere  into  the 
mid  and  upper  troposphere .  The  divergent  component  vas  the 
major  contributor,  as  the  advection  term  vas  relatively  weak 
throughout  the  budget  volume.  Processes  at  radius  10  are 
similar  vith  considerably  lover  magnitudes. 

G.  SOURCES  AND  SINKS 

Sources  cf  vcrticity  cccur  from  the  divergence  and  tilt¬ 
ing  terms,  vhile  frictional  dissipation  is  a  vorticity  sink. 
Generation  cf  vcrticity  by  horizontal  divergence  is  the 
fluid  analog  cf  the  change  in  angular  velocity  of  an  object 
when  it  changes  its  moment  of  inertia  vhile  conserving  angu¬ 
lar  momentum.  If  horizontal  divergence  occurs  and 
circulation  is  conserved,  the  area  enclosed  by  a  fluid  par¬ 
cel  vill  increase  and  the  average  vorticity  of  the  parcel 
must  decrease.  Thus,  the  pattern  of  low-level  convergence 
and  upper-level  divergence  in  a  developing  storm  means  there 
is  a  source  cf  vorticity  at  low  levels  and  a  sink  in  the 
upper  troposphere  (Fig,  42a  and  b) .  This  term  is  mathemati¬ 
cally  the  same  as  that  which  appears  in  the  lateral 
transport  partition  (se9  Section  B,  Table  V)  . 
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An  important  feature  tc  notice  is  the  relatively  large 
magnitude  of  the  maxima  in  both  regimes  compared  tc  the 
total  lateral  transport.  These  patterns  suggest  chat  the 
divergence  term  plays  the  dominant  role  in  the  budget.  In 
the  lower  troposphere,  convergence  overwhelms  all  other 
processes,  as  is  expected  from  a  scale  analysis  of  the 
general  vorticity  eguation. 

The  tilting  term  (Fig.  47a  and  b)  arises  from  vertical 
vorticity  components  generated  by  the  tilting  of  horizon¬ 
tally  oriented  vorticity  elements  by  a  non-uniform  vertical 
motion  field.  Strong  vertical  shear  and  a  horizontally  var¬ 
ying  omega  field  are  pr ereguisites  for  contributions  from 
this  source.  It  is  apparent  (from  Fig.  47a  and  b)  that  the 
products  of  vertical  shear  and  the  change  in  the  vertical 
motion  fields  in  the  horizontal  are  small  compared  to 
lateral  and  vertical  transports. 

Frictional  dissipation  (Fig.  48)  is  assumed  to  occur 
only  in  the  surface  boundary  layer  and  is  parameterized 
using  a  stability  dependent  scheme  (Johnson  and  Downey, 
1976).  The  calculated  values  were  about  1.5  times  greater 
than  those  observed  by  Ccnant  (  1982),  as  is  expected  from 
the  larger  intensity  cf  the  disturbance  studied  here.  The 
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maximum  at  1506  GOT  has  a  magnitude  very  close  to  the  value 
of  the  low-level  maximum  in  the  total  lateral  transport 
(Fig.  39a).  These  results  suggest  that  frictional  effects 
assume  a  prcminart  role  in  the  dissipation  of  vorticity  in 
the  surface  layer.  Values  in  the  outer  volume  (Fig.  48b) 
suggest  minimal  frictional  influences  on  the  vorticity 
averaged  ever  a  larger  region. 

In  summary,  the  divergence  term  was  a  major  factor  in 
the  overall  budget  as  a  source  of  cyclonic  vorticity  in  the 
lower  troposphere  and  a  sink  aloft.  The  friction  term  was 
an  important  sink  of  vorticity  in  the  surface  layer,  whereas 
the  tilting  term  had  no  significant  effect  on  the  budget. 

H.  EESIDOALS 

The  residual  in  the  vorticity  budget  contains  the 
effects  of  the  c nit  ted  processes  plus  the  accumulated  errors 
in  the  calculation  of  the  resolved  terms.  Inaccuracies  in 
the  vertical  metien  and  horizontal  wind  fields  contribute  to 
the  physical  residual  components,  whereas  spatial  and  tempo¬ 
ral  finite  differencing  and  grid-point  interpolation  inaccu¬ 
racies  are  causes  for  computational  errors.  A  positive 
vorticity  residual  indicates  an  apparent  vorticity  source  in 
which  the  ebserved  vorticity  increases  at  a  point  are  larger 
than  is  estimated  by  the  computed  terms. 
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Residuals  for  the  inner  volume  (Fig.  49a)  during  the 
development  period  indicate  generally  small  values  in  the 
middle  and  upper  troposphere,  while  large  vorticity  excesses 
are  evident  at  350  mb  (1306  GMT)  and  at  275  mb  (1406  GMT). 
Since  these  maxima  occur  during  periods  of  rapid  development 
and  maximum  vertical  motion,  it  seems  likely  that  improper 
representation  of  the  intense  mesoscale  convection  process 
is  responsible.  That  is,  inaccuracies  in  the  vertical 
motion  fields  and  therefore  the  vertical  transport  could 
lead  to  underestimation  of  the  vertical  redistribution.  The 
pattern  in  Fig.  49a  is  consistent  with  vorticiry  transport 
by  convection,  which  has  been  shown  to  be  important  in  trop¬ 
ical  circulations  (Shapiro  and  Stevens,  1982).  in  addition, 
residuals  could  also  be  expected  to  be  greater  than  normal 
in  areas  of  significant  jet  streak  activity  where  the  hori¬ 
zontal  wind  fields  are  generally  the  most  inaccurate.  The 
large  negative  residual  in  the  outer  volume  (Fig.  49b)  is  a 
likely  result  of  this  effect.  Although  computational  resi¬ 
duals  are  difficult  to  isolate  in  an  analysis  such  as  con¬ 
ducted  here,  it  appears  that  contributions  from  this  source 
are  relatively  unimportant  given  the  small  residual  values 
through  most  of  the  troposphere. 
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I.  CONCLUSIONS 

In  summary,  the  increase  of  absolute  vorticity  in  the 
budget  volume  during  the  development  stage  was  shown  to 
occur  frca  a  complicated  interaction  of  processes.  The 
period  of  surface  layer  circulation  organization  associated 
with  incipient  development  is  distinctly  different  from  the 
period  of  explosive  development  in  which  the  circulation 
increases  throughout  the  troposphere. 

During  the  incipient  stage  (period  1206-131  8),  the  lower 
troposphere  circulation  increases  ware  shown  to  occur  prima¬ 
rily  from  the  source  divergence  term,  as  the  storm  organized 
along  the  tarcclinic  zone.  Lateral  transport  increased  as 
the  circulation  developed  and  contributed  to  small  net  vcr- 
ticity  increases  since  the  frictional  dissipation  was  also 
enhanced. 

Mid-troposphere  vorticity  values  increased  in  the  earli¬ 
est  period  due  to  weak  lateral  transport  in  the  form  of  PVA 
from  the  short  wave  established  ever  Japan  (discussed  in 
detail  in  earlier  sections) .  This  mid-level  vorticity 
import  acted  to  enhance  the  low-level  organization  in  the 
earliest  stages.  At  later  times  prior  to  1318  GMT,  the 
divergence  of  the  vertical  transport  was  seen  to  account  for 
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a  significant  portion  of  the  mid-level  vorticity  increases. 
An  increase  in  the  deep  divergence  layer  above  the  LND 
served  to  offset  a  portion  of  the  mid-level  increases. 

In  the  upper  troposphere,  vorticity  values  increased  as 
the  jet  core  propagated  into  the  budget  volume.  Additional 
increases  cf  saaller  magnitudes  were  evident  from  vertical 
redistribution  from  the  surface  layer.  These  two  processes 
combined  tc  mere  than  offset  negative  vorticity  due  to  the 
increased  divergence.  Thus,  upper  level  forcing  was 
available  for  surface  development. 

The  most  significant  changes  to  account  for  the  explo¬ 
sive  development  after  1318  GMT  occurred  in  the  upper 
troposphere.  The  small  positive  lateral  transport  was  a 
balance  between  the  mean  mode  (divergence)  and  eddy  mode 
(advection)  contributions.  Large  amounts  of  cyclonic  shear 
vorticity  were  advected  into  the  budget  volume  after  1306 
GMT  as  the  jet  streak  migrated  into  the  budget  volume.  The 
corresponding  divergence  increases  occurred  12  h  later, 
which  served  tc  promote  the  explosive  low-level  development. 
The  maximum  divergence  established  a  time  period  later 
served  to  ventilate  the  system  and  induce  the  maximum  mass 
outflow  aloft,  increasing  the  surface  circulation,  and 
promoting  the  explosive  development. 
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Mid-troposphere  vorticity  increases  after  1318  GMT 
resulted  primarily  from  the  divergence  of  the  vertical 
transport  and  the  upward  development  of  the  mid-level 
diffluent  trough  from  thermal  advection.  The  satellite 
sequence  in  Chapter  4  indicates  that  convection  becomes 
intense  and  widespread  during  the  period  due  to  vertical 
stretching  in  the  mid-trcposphere  induced  by  the  increased 
upper-level  divergence.  Vertical  stretching  acts  to  cool 
the  middle  troposphere  and  destabilize  the  vertical  columns. 

In  the  lower  troposphere,  intense  circulation  increases 
resulted  frcm  the  induced  upper-level  forcing  associated 
with  the  jet  streak  movement  into  the  budget  volume.  The 
net  vorticity  increases  results  from  the  excess  of  the 
lateral  transport  over  the  frictional  dissipation.  Although 
the  frictional  effects  increased  at  about  the  same  rate  and 
to  the  same  magnitude  as  the  lateral  transport,  there  was  a 
lag  of  12  b.  When  the  upper  level  forcing  and  ventilation 
decreased  as  the  surface  center  was  removed  from  the  jet 
influence,  the  system  dissipated  after  1506  GMT. 
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VII. 


CONCLUSIONS  AND  RECOMMENDATIONS 


A  syncptic  investigation  was  conducted  of  a  western 
Pacific  case  of  ertr atropical  explosive  cyclogenssis  during 
12-17  January  1979.  Storm  related  mass  and  vorticity  budg¬ 
ets  were  computed  for  the  cyclone  using  the  2CMWF  FGGE 
analyses.  Overall  storm  development  features  were  examined 
and  discussed. 


The  following  significant  results  ware  developed: 


•  Maximum  circulation  increases  due  to  rapid  pressure 
falls  match  well  temporally  and  spatially  to  maxima  in 
the  surface  layer  inward  mass  transport  (convergence) 
and  upper-level  outward  transport  (divergence). 


•  PVA  in  the  mid- troposphere  in  the  form  of  short  wave 
interaction,  in  combination  with  thermal  advection 
effects,  enhances  the  organization  of  the  low-level 
circulation  in  the  initial  stages  of  development; 


•  Th®  initial  short  wave  interaction  decreases  after  the 
initial  period  while  surface  development  continues. 
Explosive  development  occurs  coincidentally  with  the 
migration  of  a  strong  jet  streak  into  the  budget  vol¬ 
ume,  which  provides  the  necessary  upper-level  forcing 
from  divergence  and  shear  vorticity  advection.  The 
lack  of  an  upper-level  migratory  short  wave  evident 
during  explosive  development  confirms  Pettersser.  type-A 
development; 


•  The  stability  analysis  reveals  trends  consistent  with 
numerical  simulations  (Sandgathe,  1981).  Stability 
decreases  rapidly  during  surface  development  as  cola 
air  was  drawn  over  the  warm  sea  surface.  Stability 
values  increased  after  maximum  storm  development,  as 
cold  air  advecticn  dominated  the  lower  troposphere. 


•  Mass  transport  radial  cross  sections  suggest  explosiv: 
development  characteristics  similar  to  those  observed 
in  tropical  storm  intensification  with  the  followina 
similarities:  1.  lew  static  staoility  in  the  surface 

layer.  2.  the  availability  of  an  oceanic  heat  source, 
3.  strong  outflow  (divergence)  aloft,  and  4.  strong 
convection  tear  the  inner  radii; 
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•  The  FGGE  data  set  was  found  to  be  less  than  optimum  for 
supporting  6-h  QLD  budget  analyses.  The  0600  and  1800 
GMT  time  periods  contain  spurious  potential  temperature 
biases  affecting  all  of  the  budget  calculations.  The 
problem  also  affects  both  initialized  and  unimtalized 
fields  which  sugaests  inaccurate  analysis  inputs  to  the 
model  first  guess. 


Recommendations  for  further  research  include: 


•  The  O'Brien  correction  scheme  applied  to  the  horizontal 
mass  transport  should  be  examined  to  establish  a  more 
optimum  representation  of  the  vertical  motion  fields; 


•  Further  efforts  should  be  made  to  fully  assess  the 
impact  cf  the  FGGE  0600  and  1800  GMT  data  problems  and 
their  effects  on  represented  fields  in  the  FGGE  data 
set. 


•  Energy  diagnostic  budget  calculations  should  be  com- 

fleted  on  this  cyclone  to  isolate  the  contribution  of 
hernial  advection,  sensible  heat  and  latent  heat  in  the 
explosive  development.  Correct  parameterization  of 
these  inputs  could  serve  as  a  basis  for  further  compar¬ 
isons  between  explcsive  extratropical  and  tropical 
storm  formation. 


Accurate  and  reliable  forecasts  of  maritime  explosive  cyclo¬ 
genesis  will  remain  a  difficult  problem  until  observational 
networks  are  improved,  latent  and  sensible  heat  flux  inputs 
to  storm  dynamics  are  better  parameterized,  and  vertical  and 
horizontal  model  resolutions  are  improved. 
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APPENDIX  A 


FGGE  DATA  EIAS  CHARACTERISTICS 

The  following  inforiaation  is  presented  to  document  the 
6- h  FGGE  data  problem  encountered  in  this  study.  As  indi¬ 
cated  in  Chapter  5,  a  significant  high  potential  temperature 
bias  was  observed  in  the  850-500  mb  layer  at  1800  GMT  14 
January.  Several  additional  fluctuations  were  observed  to 
occur  at  0600  and  1800  GMT  on  other  days.  The  nature  of 
these  temperature  biases  radically  affected  the  budget  cal¬ 
culations  and  precluded  the  planned  use  of  a  6-h  interval 
format.  The  nature  of  the  effects  is  shown  in  Figs.  50  and 
51.  The  top  portion  of  Fig.  50  is  a  6-h  representation  of 
the  vorticity  time  tendency  term  at  radius  six,  whereas  the 
bottom  field  is  the  inner  volume  budget  residual.  The  sig¬ 
nificant  6-h  oscillations  in  both  fields  are  indicative  of 
the  difficulty  encountered  using  QLD  budget  techniques  on 
this  particular  FGGE  data  set. 

An  even  more  striking  example  of  the  complexity  of  the 
6-h  data  problem  is  shown  in  the  inner  volume  stability 
trace  (Fig.  51).  Notice  the  oscillation  between  maxima  and 
minima  indicating  the  alternate  warming  and  cooling 
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tendencies  in  the  6-h  fialds.  It  is  obvious  that  a  synoptic 
•time  scale  interpretation  cf  this  figure  is  virtually  impos¬ 
sible.  The  1300-1500  GST  period  is  especially  interesting. 
The  0600  and  1800  GMT  stability  values  are  markedly  larger 
than  at  OOoO  and  1200  GMT.  The  increased  stabiltiy  during 
these  times  occurs  from  warm  biases  in  the  mid-level 
potential  temperatures,  as  shown  in  Fig.  27a  and  b. 

To  determine  possible  causes  for  this  6-h  data  problem, 
specific  mid-tropospheric  observational  data  available  as 
input  to  the  FGGE  analysis  scheme  at  1418  GMT  was  studied. 
The  observations  consisted  solely  of  satellite  input  (micro- 
wave  and  infrared  wavelengths).  Close  scrutiny  of  the  sat¬ 
ellite  input  in  comparison  with  the  FGGE  gridpoint  estimates 
revealed  interesting  discrepancies.  The  final  FGGE 
potential  temperature  analysis  was  5-7  K  higher  than  the 
associated  satellite  data  in  the  southwest  quadrant  of  the 
budget  volume.  In  all  other  regions  the  soundings  and  the 
FGGE  fields  matched  very  closely.  In  addition,  potential 
temperature  biases  of  smaller  magnitudes  were  evident  at  the 
1800  GMT  time  period  on  ether  days  (Fig.  27)  . 

Further  insight  is  gained  by  an  inspection  of  the  unini¬ 
tialized  height  fields  (through  the  thickness  analyses  in 
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Figs.  16c,  19c  and  22c)  .  These  fields  also  show  fluctua¬ 
tions  of  the  same  scale  occurring  at  the  1800  GMT  time 
frame,  as  this  thickness  field  departs  markedly  from  the 
preceding  ard  subsequent  analyses.  It  is  obvious  that  the 
6-h  data  problem  affects  both  initialized  (potential  temper¬ 
ature)  and  uninitialized  (height)  fields.  The  evidence  pre¬ 
sented  above  suggests  that  a  problem  existed  in  the  analyzed 
potential  temperature  fields  used  in  the  first  guess,  and 
that  the  nature  cf  the  problem  was  not  identified  during  the 
FGGE  data  post-processing  evaluations. 

Further  support  of  these  conclusions  is  presented  by 
Paegle  (1983),  who  summarized  features  of  the  FGGE  data  set 
in  tropical  regions.  Lower  tropospheric  temperature  biases 
on  the  order  cf  2-3  degrees  were  found  for  some  time  peri¬ 
ods,  which  were  presumably  from  the  lack  of  conventional 
data  input  at  the  1800  GMT  times.  In  addition,  the  ECMWF 
forecast  model  was  shewn  tc  be  insensitive  to  diurnal  temp¬ 
erature  changes  over  land.  To  what  degree  this  is  evident 
over  oceanic  areas  is  still  uncertain. 

In  summary,  it  is  believed  that  the  spurious  warming  at 
1800  GMT  between  850  and  500  mb  is  due  to  a  problem  with  the 
FGGE  analysis  scheme  rather  than  a  lack  of  quality 
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observations. 


The  iapact  of  this  assessment  on  the  results 
of  this  thesis  cannot  be  fully  ascertained  from  the  prelimi¬ 
nary  cbser vat icns  discussed  here.  Detailed  investigation 
into  specific  analysis  fields  used  in  the  first  guess  must 
be  dene  level  by  level  tc  fully  assess  the  impact  of  this 
6-h  data  problem. 
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AFPENDI X  B 


TABLES 


1 0  ib 
20  mb 
30  mb 
50  mb 
70  mb 
lOOrnt 
15  Omb 
200mt 
250mt 
30  Omt 
40  0mfc 
500mt 
70  Om  fc 
850mfc 
lOOOmb 


IABLE  I 

Availability  of  FGGE  Data 
(From  Eengsston,  1982) 
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TABLE  II 

Listing  of  Analyses  and  Prognoses  Utilized 


LEVEL 

sea-level 

sea-level 

sea-surface 

500-mb 


(FGGE*  1  , NMC=2  ,  FN0C=3) 


TYPE 

ANAL/PROG 

ORIGI 

pressure 

anal 

1,2,3 

pressure 

36-h  prog 

2,3 

te  aperature 

anal 

3 

he  ights 

anal 

1,2,3 
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TABLE  IV 

QLD  Mass  Budget  Equations  in  Isobaric  coordinates 


(After  Wash,  1978) 

The  Definition 

f  1  2 

M  m  I  —  r  sin8  dadBdp 
V 

P 

where  f  =  1  . 
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TABLE  V  (Cont.) 


Section  B 


The  partitioned  form  of  the  vorticity  budget  equation  is 


«(C) 

a 

3t 


+  LT  (C  ) 
a 


DVT  (S  ) 
.  a 


+  S  (C  ) 
a 


i 

vertical 

divergence 

- r 

vertical 

advection 

I 

mean  mode 

"  \ 

eddy  mode 

i 

horizontal 

divergence 

i 

horizontal 

advection 

divergence 

term 

tilting 

term 

The  above  partitions  make  use  of  Stokes'  theorem 


/ 


U 


m  dJl  = 


C  U  dA  , 
a  — 


and  the  division  of  total  flux  (U  £  )  into  divergent  and  advective 
*  ~  a 

components , 


v*?  u  -  5  (V*u)  +  u*v  ;  . 

a  ~  a  -  -  a 
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1EPENDI X  C 


FIGURE  S 


Figure  1.  The  Overall  Data  Management  During  FGGE  (From 
Bengsstoa,  1982). 
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Figure  2.  Stcra  Track  Positions  and  SLP  Analyses.  (A) 

Stora  Positions  in  12-h  Incraaents.  (B)  Sea 
Level  Pressure  Analyses  in  12-h  Xncreaents.  Tines 
12^0000  and  1200  refer  to  0000  SMT  12  January 
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6.  DBS E  Infrared  Satellite  Imagery  for  0041  GHT  12 
January  1979.  labels  are  (A)  Incipient  System 
(B)  Secondary  Developing  System  (C)  PVA  Area 
upstream. 
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Figure  8.  Sane  as  Figure  5  except  for  2341  GMT  12  January 
1979.  Label  &  is  Incipient  System. 
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Figure  9. 


Sage  as  Pigure  6  except  for  234  1  GST  12  January 
1979.  Label  A  is  Incipient  Systen. 
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Fscc  SST  Analysis  Valid  1200, GHT  14  Januar 
1S79.  "X"  aarks  cyclone  position  at  0000 

January.  Isotherms  in  degrees  Celsius. 


Figure  12.  Saae  as  Figure  5  except  for  1405  GMT  13  January 
1S79.  Label  A  Indicates  Jet-Frontal 
Intersection. 


ill 
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Figure  15. 


Saae  as  Figure  6  except  for  0023  GST  14  January 
1979.  Label  i  Indicates  Feeder  Bands. 
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Figure  21.  Sane  as  Figure  6  except  for  2124  GMT  14  January 

1  5  /  • 
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Figure  25.  Sane  as  Figure  6  except  for  1005  GMT  16  January 
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Figure  27. 


Area-Averaged  Potential  Temperature  Time 
Sections  (A)  Radius  6  (B)  Radius  10.  Va 
in  Degrees  Kelvin.  Time  1200  Refers  to 
12  January. 
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Figure  32.  QIC  Derived  Vertical  Velocities  (A)  Radius  6  (B) 
Radius  10.  Onits  are  ir  mb/1000-sec.  Time 
Periods  as  in  Figure  29. 
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Figure  35.  Budget  Voluae  Absolute  Vorticiry  for  Radius  6 

and  Radius  10  for  Two  Tines.  Time  1300  z  Refers 
tc  0000  GST  13  January. 
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Total,  Mean  and  Eddy  Mode  Lateral  Vorticity 
Transport  During  Maximum  Development  .  Time  a 
in  Figure  35.  units  as  in  Figure  36. 
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Figure  46 


Saae  as  Figure  37  except  for  Vertical  Vorticitv 
Advection * 
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Figure  47.  Sane  as  Figure  37  except  for  the  Tilting  Tera 
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ure  50. 


FGGS  6-h  Radius  6  Vorticity  Budget  Fields  for 
(A)  Tiae  Tendencies  (B)  Budget  Residuals. 
Period  1203  Refers  to  0000-0600  GST  12  January 
and  is  Valid  for  0300  GMT.  Contours  as  in 
Figure  37. 
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